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FOREWORD 

As  naval  operations  move  into  littoral  enviitnunents  and  address  Tactical  Ballistic  Mis¬ 
siles  Defense  (TBMD),  the  demand  for  shipboard  sensor  resoiirces  will  increase  significantly. 
FHirthermore,  threat  trends  indicate  that  the  conditions  under  which  hostile  targets  can  be 
engaged  successfully  are  beccmiing  more  difficult  to  achieve.  In  contrast  to  addressing  these 
challenging  problems  by  focusing  on  the  proctirement  of  more  information  through  more 
efficient  hardware  and  signal  processing  or  additional  sensors,  this  report  focuses  on  the  effi¬ 
cient  utiliaation  of  the  sensor  information  through  modem  tracking  algorithms  with  adaptive 
revisit  times  of  a  phased  array  radar.  This  adaptive  tracking  technology  is  presented  in  the 
form  of  a  benchmark  problem. 

While  the  problem  of  tracking  maneuvering  targets  has  been  studied  extensively,  the 
authors  were  the  first  to  define  in  the  literature  standard  problems  that  could  be  used  for 
comparison  and  evaluation  of  tracldng  algorithms.  The  lack  of  benchmark  problenu  has  hin¬ 
dered  the  progress  of  the  target  tracking  community  with  respect  to  individual  applications 
of  research  to  many  “real  world'  problems.  Thus,  many  researcbsi '  and  their  solutions  have 
not  been  considered  during  the  development  and  demonstration  of  systems  with  advanced 
technology.  In  order  to  address  this  lack  of  benchmark  problems,  the  authors  organized 
invited  sessions  at  the  1994  and  1995  American  Control  Conferences  and  invited  researchers 
to  present  solutions  to  benchmark  problems.  The  benchmark  problems  included  many  real 
world  issues  such  as  finite  sensor  resolution,  track  initiation,  beam  pointing  control,  false 
alarms,  and  Electronic  Counter  Measures  (ECM).  The  benchmark  problem  presented  in 
this  report  was  developed  from  the  lessons  learned  at  the  American  Control  Conferences. 

A  computer  program  that  will  simulate  the  adaptive  tracking  of  maneuvering  targets  as 
a  benchmark  problem  is  included  on  a  diskette  so  that  other  researchers  can  implement 
and  evaluate  their  algorithm  with  minimal  effort.  This  benchmark  has  been  accomplished 
through  fimding  from  the  NSWCDD  In-house  Laboratory  Independent  Research  (ILIR)  Pro¬ 
gram  sponsored  by  the  Office  of  Naval  Research. 

This  report  has  been  reviewed  by  Dr.  T.  R.  Rice,  Tedmicai  Lead,  Target  IVacking  and 
Signal  Processing;  and  R.  N.  Cain,  Head,  Combat  Systems  Technology  Group. 

Approved  by: 

MARY  E.  LACEY,  Head 

Systems  Research  and  Technology  Department 
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CHAPTER  1 
INTRODUCTION 


While  the  problem  of  tracking  maneuvering  targets  has  been  studied  extensively  [1-5], 
no  standard  ot  bendimark  problems  had  been  identified  in  the  litaature  for  comparison 
and  evaluatkm  of  pr^^Kwed  algorithms  until  [6,7].  While  design  objectives,  operating  condi¬ 
tions,  and  system  constraints  differ  significantly  between  tracking  problems,  clearly  defined 
benchmark  problems  in  the  areas  of  target  tracking  such  as  data  association,  multiple  target 
tracking,  and  tracking  maneuvering  targets  have  been  found  to  be  very  helpful  in  the  assess- 
jrent  and  comparison  of  existing  algorithms.  The  lack  of  benchmark  problems  has  hindered 
the  progress  of  the  target  tracking  community  with  respect  to  many  "real-world”  problems, 
and  thus,  many  researchers  and  their  sdutions  have  been  omitted  from  consideration  during 
the  devdopment  and  demonstration  of  systems  with  advanced  technology. 

The  first  benchmark  [6]  included  many  of  the  features  of  an  actual  phased  array  radar 
and  realistic  scenarios  as  well  as  many  of  the  restrictions  that  occur  in  a  "real-world”  en¬ 
vironment.  The  testbed  simulation  for  the  first  benchmaric  problem  included  the  effects  of 
target  amplitude  fluctuations,  beamshape,  missed  detections,  finite  sensor  resolution,  tar¬ 
get  maneuvers,  and  track  loss.  The  tracking  requirements  were  specified  in  terms  of  limits 
on  the  position  and  maneuverability  of  the  targets  rather  than  by  a  set  of  scenarios.  The 
"best”  tracking  algorithm  was  determined  to  be  the  one  that  minimized  the  average  number 
of  radar  dwells  while  satisfying  a  constraint  of  4  percent  on  the  maximum  number  of  lost 
tracks.  Using  a  filter  performance  criterion  based  on  system  performance  helped  to  focus 
the  efforts  of  the  researchers  onto  those  issues  that  are  important  to  the  system  designers 
and,  thus,  illustrated  the  benefits  of  modem  tracking  (i.e.,  data  processing)  through  the 
benchmark.  However,  the  bendimark  of  [6]  did  not  include  False  Alarms  (FAs)  (i.e.,  false 
detections)  or  Electronic  Counter  Measures  (ECM),  which  are  two  critical  elements  of  any 
"real-world”  tracking  problem  [8,9].  While  parameter  control  and  tracking  in  the  presence 
of  false  alarms  was  ocmsidered  fm  a  phased  array  radar  in  [10],  no  standard  problem  was 
available  to  the  authors  to  assess  the  performance  of  their  new  techniques  relative  to  other 
existing  techniques. 

This  report  extends  the  bendunark  of  [6]  to  include  the  effects  of  FAs  and  ECM.^  The 
indusion  of  FAs  is  accompanied  with  multiple  radar  waveforms  so  that  the  waveform  energy 

^  While  this  report  has  the  same  objective  as  [7],  this  report  aod  the  eonespondiof  simulation  program 
include  many  impiovemaats  over  that  of  [7]. 
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Figure  1.1  Illustration  oi  Beam  Pointing  Ccmtrol  of  a  Phased  Array  Radar 

can  be  coordinated  with  the  tracking  algorithm.  The  ECM  includes  Range  Gate  Pull  Off 
(RGPO)  on  the  target  and  a  Stand<^  Jammer  (SOJ)  broadcasting  wideband  noise.  Figure 
1.1  gives  an  illustration  ci  the  problem,  where  the  helicopter  is  a  SOJ  and  the  dotted  aircraft 
illustrates  the  effects  of  a  RGPO.  If  the  pointing  of  the  radar  beam,  denoted  by  the  strught 
solid  lines,  is  bad,  the  target  will  not  be  detected.  For  this  benchmark  the  “best”  tracking 
algorithm  is  the  one  that  minimizes  a  weighted  average  of  the  radar  energy  and  radar  time, 
while  satisfying  a  constraint  of  4%  on  the  maximum  number  of  lost  tracks. 

Figure  1.2  shows  a  block  diagram  of  the  MATLAB^  simulation  program  for  the  bench¬ 
mark,  while  additional  documentation  on  the  use  of  the  simulation  program  is  given  in 
Appendix  B.  Each  benchmark  participant  codes  his  tracking  algorithm  in  the  block  entitled 
“IVacking  Algorithm,”  which  is  given  the  range,  bearing,  and  elevation  of  the  initial  detec¬ 
tion  of  the  target.  For  each  experiment,  the  tracking  errors,  radar  energy,  and  radar  time 
are  saved.  After  the  last  experiment  of  the  Monte  Carlo  simulation,  the  average  tracking 
errors,  average  radar  energy  per  second,  and  average  radar  time  per  second^  are  computed 
for  maintained  tracks  and  the  percent  of  lost  tracks  is  also  computed.  A  constraint  of  4% 
is  to  be  imposed  on  the  number  of  lost  tracks.  A  track  is  considered  lost  if  the  distance 
between  the  true  target  position  and  the  target  position  estimate  exceeds  one  beamwidth  in 
angle  or  1.5  range  gates.  When  FAs  and  ECM  are  present  in  an  actual  radar  system,  algo¬ 
rithms  for  reacquiring  the  target  and  “coasting”  the  target  tracks  through  jamming  signals 
required  in  order  to  maintain  a  track.  However,  in  this  benchmark  problem,  the  SOJ 
power  is  limited  so  that  the  it  can  be  defeated  with  one  of  the  higher  energy  waveforms.^ 

^  MATLAB  is  a  trademark  of  The  MathWorks,  Inc.,  Natick,  MA. 

^  This  is  inversely  proportional  to  the  average  sample  period.  More  details  on  this  are  given  in  Chapter  6. 

^  This  restriction  is  made  to  limit  the  scope  of  the  problem  and  reduce  the  problems  associated  with 
declaring  a  track  as  lost. 


1-2 


NSWCDD/TR-96/10 


Since  track  initiation  is  very  expensive  with  respect  to  radar  time  and  energy,  the  emphasis 
of  this  benchmark  is  on  muntaining  triu^ks  with  an  allowance  for  track  reacquisition  within 
a  few  radar  dwells.’ 

The  radar  model  simulation  includes  the  effects  of  target  amplitude  fluctuations, 
beamshape,  missed  detections,  finite  resolution,  FAs,  an  SOJ,  and  RGPO.  The  radar  per¬ 
forms  search  dwells,  monopulse  track  dwells,  and  monopulse  pusive  dwells.  Eight  radar 
waveforms  that  differ  primarily  in  the  pulse  length  are  available  for  control  (i.e.,  selection) 
&om  the  tracking  (and  radar  management)  algorithm.  Since  a  waveform  with  a  longer  pulse 
length  provides  a  higher  Signal-to-Noise  Ratio  (SNR)  at  the  cost  of  more  radar  energy,  the 
proper  coordination  of  the  waveform  selection  with  the  tracking  algorithm  is  an  important 
“real-world”  issue  to  be  addressed  by  each  participant.  For  example,  a  higher  target  SNR 
results  in  fewer  FAs  at  a  cost  of  more  radar  energy.  RGPO  and  standoff  januning  in  the 
mainlobe  tmd  sidelobes  are  included  in  the  benchmark  problem.  In  RGPO,  the  target  under 
track  repeats  with  delay  the  radar  pulse  to  pull  the  radar  range  gate  off  the  target  so  that  no 
detection  of  the  true  target  will  occur.  The  time  delay  is  controlled  so  that  the  false  target 
is  separated  from  the  target  with  linear  or  quadratic  motion.  For  tracking  an  SOJ,  passive 
measurements  of  the  janruner  position  (i.e.,  angle  only)  are  provided  by  the  radar  at  the 
request  of  the  tracking  algorithm,  with  bearing  and  elevation  angles  for  pointing.  Thus,  the 
tracking  algorithm  will  be  required  to  track  one  target  and  one  jammer.  When  the  jammer 
is  in  the  mainlobe  of  the  radar  beam  pattern,  the  target  return  will  be  corrupted  or  hidden 
by  the  jammer  signal.  When  the  jammer  is  in  one  of  the  sidelobes,  the  SNR  for  the  target 
will  be  reduced.  The  initial  bearing  and  elevation  of  the  jammer  will  be  given  to  the  tracking 
algorithm. 

The  targets  exhibit  Radar  Cross  Section  (RCS)  fluctuations  according  to  the  Swerling  3 
type  and  perform  as  much  as  7  g  of  lateral  acceleration  and  2  g  of  longitudinal  acceleration. 
Target  range  can  vary  from  20  to  100  km,  while  the  target  elevation  angle  can  vary  from  2° 
to  30**.  Since  only  one  face  of  the  phased  array  radar  is  used,  the  bearing  (or  azimuth)  of  the 
target  will  be  confined  to  ±60**.  The  average  RCSs  of  the  targets  are  large  enough  so  that 
an  average  SNR  of  18  dB  is  achievable  with  the  highest  energy  waveform.  The  SOJ  remains 
at  ranges  near  150  km  and  performs  less  than  2  g  of  acceleration. 


Results  of  Previous  Benchmark  Efforts 

An  invited  session  was  organized  at  the  1994  American  Control  Conference  (ACC)  with 
the  first  benchmark  as  the  theme  for  the  session.  The  benchmark  problem  was  presented, 
along  with  four  different  approaches  to  the  problem.  The  results  of  the  invited  session  [11-13] 
and  other  studies  such  as  (14)  provided  relative  comparisons  of  different  tracking  algorithms 
on  the  first  benchmark.  The  a-  P  filter  provided  an  average  sample  period  of  0.85  s,  while 


’  Reacquisition  is  allowed  if  it  can  be  acconiplished  before  the  estimation  errors  exceed  the  criteria  for  a 
track  being  considered  lost. 
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a  staadard  Kalman  filter  with  a  nearly  constant  velocity  motion  model  provided  an  average 
sample  period  of  about  1  s.  Thus,  for  an  order  of  magnitude  increase  in  computations, 
the  Kalman  filter  provided  only  an  20%  increase  in  the  average  sample  period.  A  two- 
model  Interacting  Multiple  Model  (IMM)  algorithm  provided  an  average  sample  period  of 
1.3  8,  while  a  three-model  IMM  algorithm  provided  an  average  sample  period  of  1.5  s. 
Incorporating  adaptive  revisit  times  with  a  three-model  IMM  algorithm  provided  an  average 
sample  period  of  2.3  s,  for  about  1.6  orders  of  magnitude  increase  in  computations  relative  to 
the  a  —  P  filter.  An  H- Infinity  filter  solution  was  originally  included  in  the  invited  session  at 
the  1994  ACC,  but  the  authors  withdrew  the  paper  when  they  concluded  that  the  H-Infinity 
filter  provided  no  significant  advantage  over  the  Kalman  filte.  [15]. 

The  panel  discussion  held  during  the  invited  session  at  the  1994  ACC  revealed  the  im¬ 
portance  of  tracking  maneuvering  targets  in  the  presence  of  ECM.  Thus,  a  second  invited 
session  was  orgemized  for  the  1995  ACC  with  a  new  benchmark  [7]  that  included  fr  'se  alarms 
and  ECM.  However,  the  difficulty  of  the  new  benchmark,  discrepancies  in  the  calculation 
of  the  radar  energy,  and  deficiencies  in  the  simulation  program  resulted  in  no  quantitative 
results  for  comparison  of  the  performances  of  the  different  algorithms.  Qualitatively,  the 
IMM  algorithm  and  Multiple  Hypotheses  Tracking  (MHT)  yielded  comparable  results,  with 
the  MHT  algorithm  being  1  to  2  orders  of  magnitude  costlier  in  computations. 

Many  of  the  deficiencies  associated  with  the  benchmark  problem  at  the  1995  ACC  were 
corrected  to  produce  the  benchmark  problem  described  in  this  report.  Reference  [16]  presents 
an  IMMPDAF  solution  to  this  benchmark  problem,  while  [17]  presents  an  adaptive  Kalman 
filter  solution.  The  adaptive  Kalman  filter  provided  an  average  sample  period  near  1  s,  while 
the  IMMPDAF  provided  an  average  sample  period  near  2  s.  The  adaptive  Kalman  filter 
also  required  about  30%  more  energy  than  the  IMMPDAF.  To  date,  no  results  for  MHT  are 
available  to  assess  its  performance  relative  to  these  two  solutions. 


Orgsuiiasation  of  the  Report 

This  report  is  organized  as  follows.  Chapter  2  describes  the  radar  model,  while  Chapter  3 
discusses  the  SOJ  and  RGPO.  The  target  trajectories  are  presented  in  Chapter  4,  along 
with  the  average  RCS  for  each  target  and  the  timing  for  the  use  of  the  RGPO.  Chapter  5 
discusses  the  inputs  and  outputs  of  the  tracking  algorithm.  The  criteria  for  evaluating  the 
performance  of  the  tracking  algorithm  are  given  in  Chapter  6,  and  concluding  remarks  are 
given  in  Chapter  7. 
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CHAPTER  a 
RADAR  MODEL 


The  rader  is  a  4  GHc  phased  array  using  amplitude-comparison  monopulse  with  uniform 
illumination  across  the  arr^.  Each  radar  dwell  consists  of  one  phaae/frequency  discrete- 
coded  pulse  [18].  The  range  gate  is  approximately  1575  meters  for  a  track  dwell  and  10  km 
for  a  search  dwell,  wi'.h  the  number  of  range  bins  varying  from  70  to  444  depending  upon  the 
waveftmn  and  dwell  ^ype.  The  radar  beam  is  quasi- circular  with  the  beamwidth  increasing 
as  the  beam  is  steered  off  the  broadside  direction.  The  radar  beam  has  a  3  dB  beamwidth 
of  —  2.4^  on  broadside  (i.e.,  normal  to  the  face  of  the  array)  and  Bgyy  4-5**  at  a 
broadside  angle  of  60".  The  two-way  radar  beam  has  a  3  dB  beamwidth  of  Bgw  — 
on  broadside  and  0bw  ~  3*2°  at  a  broadside  angle  of  60".  The  beam  is  pointed  to  the 
commanded  direction  at  the  next  available  sample  period  since  the  target  trajectories  are 
stored  in  the  data  file  at  20  Hx.  The  minimum  time  period  between  sets  of  radar  dwells  is 
restricted  to  0.1  s  (i.e.,  10  Hz).  A  set  of  radar  dwells  may  consist  of  as  many  as  five  dwells 
that  are  requested  simultaneously  by  the  tracking  algorithm.  Each  radar  dwell  requires  0.001 
s  of  radar  time  to  be  accomplished.  However,  for  the  purposes  of  simulation,  all  of  the  dwells 
in  a  set  are  modeled  as  occurring  at  the  same  instant  of  time  with  respect  to  the  target  state. 
For  a  given  detection  threshold  selected  by  the  participant's  radar  management  algorithm, 
the  radar  mode!  will  return,  for  each  detection  t  of  each  dwell  of  the  set,  the  following:  the 
SNR  bearing  monopulae  ratio  r^,  elevation  monopulse  ratio  r]^,  and  range  r][,  where 
k  denotes  the  dwell  time,  which  is  assumed  to  be  equal  for  all  dwells  in  the  set.  The  sum 
and  difference  channels  used  to  produce  the  angular  measurements  are  corrupted  with  white 
Gaussian  errors.  Note  that  the  errors  in  the  monopulae  ratio  ate  not  purely  Gaussian. 

Baric  Radar  Equation 

The  basic  radar  equation  for  a  single  pulse  and  target  at  range  R  is  given  in  [19]  as 

„  P,G,GA‘F?F;<rt  fG^(R)\ 

- ■  ■  l-ftT-j  (21) 
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where 


Pr  = 
Pt^ 

Gr  = 
A  = 
Ft  = 

<7o  = 

i-M  = 

G»u{R)  = 


Received  power  from  the  t&rget 
Trursmitter  power 

Transmitter  antenna  gain  in  the  direction  of  the  target 

Receiver  antenna  gain  in  the  direction  of  the  target 

Transmit  wavelength 

Transmitter  propagation  factor 

Receiver  propagation  factor 

Target  RCS 

Total  losses  for  the  radar  system 
Sensitivity  Time  Control  (STC)  gain 


STC  is  employed  to  prevent  saturation  of  the  receiver  by  close-in  scattering  and  to  prevent 
short-range  detection  of  objects  with  a  small  RCS.  For  this  benchmark,  the  STC  gain  is 


defined  as 


G,u{R)  = 


R  ^  R»ie 


R  <  R*u 


(2.2) 


where  R,tc  ~  30  km.  The  antenna  gain  along  the  boresight  of  the  beam  is  computed  in  terms 
of  the  munber  of  elements  assuming  half-wavelength  separation  distance  between  elements. 
Note  that  the  gain  becomes  wavelength  independent  in  this  approximation.  The  effective 
aperture  is  given  by 


Ae  —  £aA 


(2.3) 


where  i>  the  aperture  efficiency,  and  A  =  LgL^  is  the  actual  area  for  array  dimensions 
Lg  and  Lp.  The  gain  is  given  by 


C?  = 


4irAf 

irsA 

~xr 


-^L  L 


=  rcA^tN, 


(2.4) 


where  N,  and  Nf  are  the  number  of  array  elements  in  a  particular  dimension.  In  this 
benchmark,  55  elements  per  array  side  is  assumed  along  with  an  aperture  efficiency  of  0.5. 
Thus,  Gt  =  Gt  =  36.8  dB  for  perfect  beam  pointing. 


Pn  ~ 

2-2 


The  noise  power  is  given  by 
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where 


ibo  =  Boltzmann  constant 
To  =  Reference  temperature 
Bn  =  Noise  bandwidth 
Fn  =  Receiver  noise  figure 

The  uncompressed  (i.e.,  no  phase/frequency  coding)  SNR  K*  is  given  by 

PtGtGrX^F^F^co  (  1  \  (G^c{R)\ 


(  1  \  (G^{1 

\koToBnFn)\  B* 


V 


(2.6) 


For  a  radar  with  a  given  Ft,  the  SNR  at  the  output  of  the  matched  filter  can  be  increased 
by  extending  the  length  of  the  radar  pulse.  In  order  to  achieve  good  range  resolution  with  a 
long  pulse,  pulse  compression  is  used.  In  discrete-coded  pulse  compression,  the  pulse  width, 
Tc,  is  composed  of  Ns  subpulses  of  width  (i.e.,  Tt  =  Nsts),  where  each  subpulse  is  coded 
with  frequency  and/or  phase.  For  this  benchmark  problem,  the  radar  waveform  is  biphase 
coded  [18].  The  subpulse  width  ts  is  assumed  to  be  equal  to  the  compressed  pulse  width 
Te  and,  thus,  rj’  defines  the  recdver  bandwidth  of  the  compressed  pulse  Bn.  Modeling  the 
compressed  pulse  as  coherently  integrated  Ns  subpulses,  the  output  SNR  for  the  compressed 
pulse  is  given  by 


FGtGrX^F^F^ao  (  N,t,  \  (G,^,(R)\ 


(2.7) 


Assuming  an  RCS  of  one  square  meter  and  the  radar  parameters  of  Table  2.1,  the  radar 
equation  can  be  simplified  to 

=  (2.8) 

where  Do  *  257.6  dB  .  For  each  of  the  eight  waveforms  and  a  one-square-meter  target,  the 
SNRs  at  100  km  were  calculated  using  (2.8).  The  results  are  summarized  in  Table  2.2,  where 
Ar  K  O.bcrs  is  the  resolution  of  each  pulse,  with  c  being  the  speed  of  light. 


The  RCS  of  the  target  will  be  modeled  as  Swerling  3  type,  where  the  density  function  for 
the  RCS  is  given  by 

A<r 

/(•r)  = 


'4V€ 


2a 

^mve 


(2.9) 


with  the  average  RCS  (i.e.,  Onvt)  varying  between  target  scenarios.  The  cumulative  distri- 
butim  function  of  the  radar  cross  section  is  given  by 

2^0 


F{ao)  =  P{0<<y<<7o}  =  l-  fl-1-  exp  — 


'ave 


(2.10) 
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Table  2.1  Radw  Parameters 


lUClJ] 

imxu] 


Table  2.2  Waveform  Parameters 


Since  F((Tq)  is  uniformly  distributed  between  0  and  1,  a  sample  or  observation  of  the  RCS 
ran  be  computed  by  solving 


V  I  ‘^•wj 


where  x  is  a  random  number  that  is  uniformly  distributed  between  0  and  1. 
The  power  received  from  the  SOJ  is  given  by 


where 


Jammer  power  at  the  radar 

Radar  receiver  antenna  gun  in  the  direction  of  jammer 
Effective  bandwidth  of  the  radar  receive  filters 
Effective  bandwidth  of  the  jammer  emission 


(2.11) 


(2.12) 


NSWCDD/TR-96/10 


For  this  benchmark,  the  SOJ  power  Pj  will  be  a  multiplicative  factor  of  ‘.he  noise  power  with 
Bj  >  Bf  Thus, 

P,  =  70  GrjGMRctr)  (2.13) 

where  70  s  8  will  be  fixed  during  the  benchmark,  is  the  STC  gain  evaluated  at 

the  center  range  bin  of  the  active  dwells  and  at  unity  on  passive  dwells. 

Range  Mcaeuremente 

The  range  is  computed  as  the  range  bin  value  that  is  closest  to  the  true  range  in  order 
to  eliminate  the  need  for  processing  adjacent  bins  for  detection.  Also,  since  the  targets  are 
assruned  to  be  points  in  space  and  no  interpolation  between  bins  is  performed,  the  discrete 
coding  results  in  range  measurements  with  errors  that  are  uniformly  distributed  and  defined 
by  the  bin  resolution  (i.e.,  independent  of  SNR). 

The  center  of  the  range  gate  is  placed  at  the  predicted  range  ft  provided  by  the  tracking 
algorithm.  The  range  measurement  is  computed  as  follows.  The  prediction  error  in  range  is 
given  by 

(2.14) 

where  /Z^  is  the  true  range  of  the  target  at  time  k.  If  |ft|  >  775  m,  then  a  missed  detection 
occurs  and  none  of  the  measurements  that  are  returned  will  include  the  target.  If  the  target 
is  in  the  range  gate,  the  range  measurement  is  given  by 

rj  *  Arrd[^] +f*  (2.15) 

where  Ar  is  specified  in  Table  2.2  for  each  waveform  and  rd[  ]  denotes  the  nearest  integer. 
The  range  of  an  FA  is  assigned  to  that  range  of  the  bin  in  which  it  occurs. 

Atttanna  Gain  Patterns 

The  array  consists  of  3025  individual  elements  (i.e.,  55  elements  in  elevation  and  bearing) 
with  cosine  illumination  and  half-power  beamwidth  of  140*.  The  broadside  of  the  array  is 
directed  at  0*  bearing  and  15*  elevation.  The  normalized  voltage  patterns  (20]  are  defined  by 
voltage  pattern  in  elevation  pointed  at  and  target  at  and  H(^i«^*)i 
the  voltage  pattern  in  bearing  pointed  at  ih  and  target  at  Bl,  and  are  pven  by 

e  sin(iV52)  [siu(6|  -f  0.5ir)  sin(5i  -  O.Oir)! 

■  4^.43)  Arsio(63)  1  ij  -f-  0.55r  5,  -  0,5x  J 


(2.16) 
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where 


ir  sin(A^a3) 

’8in(ai  +  0.5ir)  _  8in(ai  -  O.Sjt)] 

(0  171 

4(1.43)  Nnn{a2) 

ai  +  0.5ir  *  oj  —  0.5ir  J 

{£.11) 

»  0.25  sin(^) 

(2.18) 

“a 

=  0.5w  sin(B|)  -  sin(5fc)j 

(2.19) 

5i 

=  0.25  ain(£{  - 

-«*) 

(2.20) 

5a 

s  0.5ir  [sin(£][ 

-  ej)  -  8in(gfc  -  e*)] 

(2.21) 

N 

=  55 

(2.22) 

e» 

=  15® 

(2.23) 

The  factor  of  1.43  has  been  included  to  give  a  normalized  sum  channel  voltage  at  boresight 
for  a  broadside  angle  of  0^. 


Channel  Voltages 

The  received  voltage  will  be  related  to  the  received  power  as  V,  s  v/27^.  In  search  or  track 
mode,  the  sum  voltage  at  time  k,  is  composed  of  the  in-phase  and  quadrature  portions. 
The  in-phase  and  quadrature  components  of  the  sum  voltage  s'^  normalized  by  the  receiver 
noise  is  given  by 


where 


•'ik  cos^i  +  risi  cos^*  -I-  V(0, 1) 

>Qk  -rl(Ei)*  sin^  -I-  rjEj  sin^^j  +  N{0, 1) 

ri  -  [^3(0. 1)  -I-  iVa(0, 1)] 

”  V  19®-2  <*B 

El  -  + *!(-«:,.  *i,) + ♦!(«,,-<) + tK-c,,  -ei,) 


K{Eiii  +  xmBiX  +  v) 


(2.24) 

(2.25) 


(2.26) 

(2.27) 

(2.28) 

(2.29) 

(2.30) 

(2.31) 

(2.32) 
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(2.33) 

(2.34) 


^*io/<»3(ejb  -  et) 

*  9,fo/cos0>k) 

Tt  =  Pulse  length 

^1^  s  Phase  of  the  voltage  of  the  target  return 
^  =  Phase  of  the  voltage  of  the  SOJ 
BI  —  Bearing  of  the  SOJ  at  time  k 
El  =  Elevation  of  the  SOJ  at  time  k 

Both  and  ^  are  independent  and  uniformly  distributed  between  0  and  2ir.  The 
and  0^  denote  the  squint  angles*  in  elevation  and  bearing,  respectively,  where 
=  0.9**  =  15.7  mrad.  Also,  N{0, 1)  denotes  a  Gaussian  random  variable  with  zero  mean 
and  unity  variance.  Let 

(*i)’ ■=(»/.)’ +  (•!}*)'  (2'») 

Then  the  observed  SNR  is  computed  by  9tk  =  0.5(sj^)^. 

Eqs.  (2.24)  through  (2.34)  are  used  to  compute  a  measurement  for  each  range  bin  by 
setting  3B  0  for  the  bins  that  do  not  include  the  target.  The  passive  measurements  are 
obtained  by  setting  Fj^  -  0  and  G,tc  **  1  for  all  range  bins.  The  observed  SNR  Rj  and  range 
are  returned  for  every  range  bin  i  with  Ri  >  RJ*,  where  Rj|*  is  specified  by  the  tracking 
algorithm  for  the  dwell  set  at  time  k.  In  other  words,  a  measurement  is  returned  for  range 
bin  t  if  R^  >  R|*.  Thus,  the  probability  of  an  FA  P/«  is  specified  by  the  tracking  algorithm. 
The  amplitudes  of  the  FAs  are  Rayleigh  distributed  with  parameter  7o(Sj[)^G(icvA)  +  1* 
Since  the  passive  measurements  are  obtained  by  setting  F^  v  0  and  Ggtc  1  for  all  range 
bins,  the  amplitude  of  the  measurements  of  a  pass  "  dwell  will  be  Rayleigh  distributed  with 
parameter  7o(Si)^  +  !• 

Figure  2.1  shows  the  one-way  beamwidth  ot  a  sum  channel  versus  the  off-broadside 
angle,  which  is  the  angle  between  the  boresigbt  of  the  beam  and  a  vector  orthogonal  (i.e., 
normal)  to  the  face  of  the  array.  The  broadside  angle  is  treated  independently  in  elevation 
and  bearing.  The  one-way  beamwidth  is  the  angle  between  the  two  half-power  points  on 
the  sum  channel  beam.  The  one-way  beamwidth  varies  from  2.4*  at  broadside  to  4.5*  at  a 
broadside  angle  of  6(F.  The  two-way  beamwidth  varies  from  1.6*  at  broadside  to  3.2*  at  a 
broadside  angle  of  60*.  Figure  2.2  shows  the  normalized  antenna  gain  versus  the  off- broadside 

*  Note  that  the  squint  eagles  an  usually  fixed  with  nspeet  to  the  olf-broadaide  angle.  The  iquint  angles 
have  been  varied  with  the  HT-broadside  angle  to  increase  the  minimum  monopulse  error  elope. 
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angle,  where  the  gain  has  been  normalized  by  the  gain  at  broadside.  Thus,  antenna  gains 
Gf  and  Gr  Me  decreased  by  0.4  dB  at  a  broadside  angle  of  60°. 

The  difference  voltages  in  elevation  and  bearing  <4  in-phase  and  quadrature 
componentb.  The  in-phase  and  quadrature  components  of  and  <lj[  are  written  as 

fn  =  rlEi  «*„) + -d},)  -  «L(-d:„d:,)  -  «!(-«;,,  -<))  co.di 

+  r{[»{(»5t.»*M) + -*!i)  -  -  *{(-«!.  -*!i)) 

+  JV(0,1)  (2.36) 

di,,  =  rSEi  [♦S(«:,.<) + •}(«,,-»!,)  -  ♦»(-*i.*i,)  -  »i"dt 

+  rj [♦{(«..<)  +  ♦’.(*'«. -fit)  -  ♦’.(-«'«■<)  -  «’»(-«..  -<)] 

+  Ar(0,l)  (2.37) 

d}, = riES  [♦K«:„<) + ♦i(-«,.<<,)  -  -fit)  -  «*(-*:..  -<)]  wi 

+ r{  [«{(<^.<!,) + ♦{(-*‘.,.«i,)  -  ♦’»(«,.  -*!,)  -  »{(-«,. 

+  f/(0,l}  (2.38) 

4»  =  rSEi  + «S(-d^,d!,)  -  ♦Ud'.,,  -dj,)  -  ♦i(.-d;,,  -d*,)] 

+  ri  [♦'»(«..<)  +  *i(-Vffit)  -  ♦S(C  -*Jf)  -  «{(-*'«.  -<)) 

+  JV(0.1)  (2.39) 

where  (he  received  errors  in  the  difference  channels  have  been  modeled  as  independent  of  the 
receiver  errors  in  the  sum  channel.  Passive  measurements  are  generated  by  setting  fj^  »  0. 

Monopulee  ProceMing 

Monopulse  processing  is  a  simultaneous  lobing  technique  for  determining  the  angular 
location  of  a  source  of  radiation  or  of  a  “target”  that  reflects  part  of  the  energy  incident 
upon  it  [2]].  The  raonopulse  ratio  is  formed  by  the  (complex)  division  of  the  difference- 
channel  voltage  phasor  by  the  sum-channel  voltage  phasor.  Monopulse  ratios  are  formed 
separately  for  the  elevation  and  bearing  directions  of  arrival.  For  every  active  or  passive 
track  dwell,  monopulse  ratios  are  returned  for  every  measurement  with  an  amplitude  that 
exceeds  the  specified  detection  threshold  9^*.  The  monopulse  ratios  are  then  used  to  estimate 
the  monopulse  error  function  value  ij  for  determining  the  measured  direction  of  arrival. 
The  processing  of  monopulse  ratios  will  be  considered  as  in  Appendix  A  for  (i)  target  only 
measurements,  (ii)  SOJ  only,  and  (iii)  target  in  the  presence  of  the  SOJ  within  the  mainlobe 
of  the  antenna  pattern. 
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The  monopulse  error  function  et  broadside  is  shown  in  Figure  2.3  along  with  the  commonly 
used  linear  approximation  to  the  monopulse  error  function,  which  is  given  by 

V »  fcm®)  -^BW  <9  <  Bbw  (2.40) 

where  is  the  average  monopulse  error  slope  for  -0.5  9bw  <  9  <  0.5  $bw-  Figure  2.4 
shows  the  monopulse  error  slopes  versus  the  off-broadside  angle.  Since  the  squint  angles  are 
increased  with  the  broadside  angle,  the  monopulse  error  slope  achieves  a  minimum  of 
24  rad'*  rather  than  12  rad~*  if  the  squint  angles  had  been  fixed  at  their  boresight  values. 
Since  the  standard  deviations  of  the  angle  measurements  are  inversely  proportional  to  km, 
maintaining  a  higher  km  provides  more  accurate  measurements  at  the  off-broadside  angles. 


The  monopulse  ratios  in  elevation  and  bearing  are  computed  as 


e  _  *'tk^Ik  *Qk*^Qk 

■  is'n)^  +  (-'<,*)» 

I  +  *Qk^Qk 


(2.41) 

(2.42) 


When  detecting  a  single  target,  the  Directions-Of-Arrival  (DOA)  estimates  in  elevation 
and  bearing  i;[*  are  given  by 


»*  >  3  dB 

qj*  =  (l  +  Rt>3dB 


(2.43) 

(2.44) 


Using  the  linear  approximation  for  the  monopulse  error  function  gives  the  measurements  of 
elevation  and  bearing  as 


4  = 


(2.45) 

(2.46) 


where 


kfn  *  48c<M(e*  -  e») 
=  48  cos(iA) 


volt 

volt  •  rad 
volt 

volt  •  rad 
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and  eit  and  bh  are  the  elevation  and  bearing  pointing  commands  from  the  tracking  algorithm. 
Estimates  for  the  standard  deviations  of  the  measurements  are  given  by 


.g  1 


tftk  >Zdh 
Rs  >  3  dB 


(2.49) 

(2.50) 


Note  that  (2.49)  and  (2.50)  differ  from  the  expression  pven  in  [21]  for  the  standard  deviation. 
The  variance  expression  in  [21]  is  based  on  an  assumption  of  {i^l  <  1  and  S2(  >  12  dB,  while 
(2.49)  and  (2.50)  are  valid  for  ji;]  <  2  and  Rs  >  3  dB.  Figure  2.5  pves  the  standard  deviation 
for  various  off-broadside  steering  angles  and  R*  =  15  and  20  dB  with  perfect  pointing  (i.e., 
dj[  s  0).  Note  that  Rs  can  vary  rather  significantly  between  consecutive  measurements  since 
the  targets  have  a  RCS  of  Swerling  3  type. 


When  processing  the  passive  measurements  of  the  SOJ,  the  Jammer-to-Noise  Ratio  (JNR) 
Jk  can  be  estimated  from  N  detections  within  a  dwell  according  to 

(2M) 

iml 
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where  is  the  observed  SNR  of  the  range  hin 

elevation  and  bearing  of  the  SOJ  are  given  by 

i  at  time  k. 

The  DOA  estimates  for  the 

i  +  y 

L  Jki 

**  * 

J  >3dB 

(2.52) 

hs] 

Hk  » 

J*  >  3  dB 

(2.53) 

where 

fi‘= 

EiLi*: 

(2.54) 

Sir? 

(2.53) 

1>« 

with  r”  and  denoting  the  monopulse  ratio  for  range  bin  t  at  time  k  for  the  elevation 
and  bearing,  respectively.  The  and  are  the  minimum  variance  estimates  of  the  means 
of  the  monopulse  ratios  for  elevation  and  bearing,  respectively.  Using  (2.52)  and  (2.53)  in 
(2.45)  and  (2.46)  provides  the  measured  elevation  and  bearing  of  the  SOJ.  Estimates  for  the 
standard  deviations  of  the  SOJ  measurements  are  given  by 


%>3dB  (2.56) 
Ji>3dB  (2.57) 


When  the  target  and  the  SOJ  are  in  the  mainlobe  of  the  antenna  pattern,  the  angular 
location  of  the  SOJ  can  be  observed  as  discussed  in  (2.51)  through  (2.57)  with  a  passive 
dwell  or  an  active  dwell,  where  only  the  range  bins  that  do  not  include  the  target  are  to 
be  considered.  Thus,  the  DO  A  estimates  for  the  target  in  the  presence  of  the  SOJ  in  the 
mainlobe  are  given  by 


> 

[  2**  J 

ri- 

•  < 

Jjk  +  1 

2«* 

fje 

fftk  >2Jk  +  2 

(2.58) 

[  ^  2»*  J 

ri- 

■js+r 

2Jtk 

Hi  » 

R*  >  27»  +  2 

(2.59) 

where  and  are  the  minimum  variance  estimates  of  the  means  of  the  monopulse  ratios 
that  do  not  include  the  target.  Eqs.  (2.58)  and  (2.59)  indicate  that  if  the  observed  SNR  is 
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12  dB  greater  than  the  JNR  the  effects  of  the  SOJ  can  be  ignored.  Note  that  should  be 
estimated  on  the  active  dwell  since  the  SOJ  power  can  fluctuate  between  consecutive  radar 
I  dwells.  Using  (2.58)  and  (2.59)  in  (2.45)  and  (2.46)  provides  the  measured  elevation  and 

bearing  of  the  target  in  the  presence  of  the  jammer.  Estimates  of  the  standard  deviations  of 
the  measurements  are  given  by 


I 


ll>e  1 

'  .  :/k+i' 

* 

.ii>*  _  ' 

7k + 1 

2»* 

where  >  2Jk  +  2.  Eqs.  (2.58)  and  2.59)  provide  estimates  of  the  DOAs  of  a  target  in  the 
presence  of  an  SOJ  and  (2.60)  and  (2.61)  provide  estimates  of  the  vsiriance  associated  with 
the  DOA  estimates.  However,  note  that  t^y  and  J  have  been  used  as  known  quantities.  Then 
a  sufficiently  large  number  of  range  bins  should  be  used  to  estimate  J  and  t)j.  A  reasonable 

A 

criteria  for  M  is  MJ  >  18  dB.  For  example,  an  SOJ  with  J  =  6  dB  would  require  about  16 
range  bins. 


I 


I 
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CHAPTER  S 

ECM  TECHNIQUES  AND  MODELING 

The  ECM  for  the  benchmark  involves  the  SOJ  discussed  in  Chapter  2  and  RGPO.  The 
times  at  which  the  SOJ  and  RGPO  impact  the  trajectories  are  discussed  in  the  next  section. 
The  SOJ  transmits  broadband  noise  toward  the  radar  and  has  the  effect  of  increasing  the 
level  of  noise.  However,  since  the  SOJ  noise  enters  the  sum  and  difference  channels  with 
the  same  phase,  the  DOA  of  the  SOJ  can  be  estimated  with  the  monopulse  ratios,  where 
the  target  amplitude  is  Rayleigh,  since  the  SOJ  noise  is  modeled  as  Gaussian.  The  SOJ 
flies  an  oval  (race  course)  holding  pattern  in  a  clockwise  direction  at  an  altitude  of  3.05  km 
and  speed  of  168  m/s.  The  SOJ  is  approximately  150  km  from  the  radar,  and  its  trajectory 
is  shown  in  Figure  3.1.  The  aircraft  flies  straight  and  level  for  a  period  of  approximately 
40  s,  upon  which  a  1.0-g  turn  is  performed  through  a  ISO-degree  heading  change.  Straight 
and  level,  nonaccelerating  flight  is  continued  for  another  40  s,  where  a  second  1.0*g  turn  is 
performed  until  the  aircraft  returns  to  its  original  position.  Then  the  cycle  starts  again  with 
another  40  s  of  level  flight  maintained.  The  SOJ  transmits  broadband  noise  that  impacts  the 
radar  with  power  no  more  than  70  times  the  receiver  noise  power,  as  discussed  in  Chapter  2. 
For  this  benchmark,  70  s  8.  Thus,  the  SOJ  will  not  completely  bide  a  target  and  it  can  be 
defeated  with  a  higher  energy  waveform  in  this  benchmark.  While  the  SOJ  energy  enters 
into  every  radar  dwell,  the  effects  of  the  SOJ  are,  for  the  most  part,  negligible  when  it  is  not 
in  the  mainlobe  of  the  antenna  pattern.  Note  that  this  is  not  the  case  for  SOJs  transmitting 
a  higher  level  of  power. 

In  RGPO,  the  target  under  track  repeats,  with  delay  and  amplification,  the  radar  pulse 
so  as  to  pull  the  radar  range  gate  off  the  target.  The  time  delay  is  controlled  so  that  the 
false  target  is  separated  from  the  target  with  linear  or  quadratic  motion.  For  the  linear  case, 
the  range  of  the  false  target  with  respect  to  the  radar  is  expressed  as 

=  +  (3.1) 

where  is  the  range  of  the  actual  target,  is  the  rate  of  pull-off,  4  >>  the  time  of  the 
radar  dwell,  and  to  is  the  initial  reference  time  of  the  RGPO  false  target.  For  the  quadratic 
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case,  the  range  of  the  false  target  with  respect  to  the  radar  is  written  as 

=  /?*  +  ^app{tif  —  to)*  (3.2) 

where  is  the  acceleration  of  pull-off  of  the  false  target.  The  measurements  of  the  RGPO 
will  be  generated  as  those  for  the  target  with  (2.26)  modified  according  to 

T«7l<^wGrtc(/**)  (3.3) 

where  71  is  the  amplification  factor  of  the  RGPO.  However,  when  the  times  of  arrival  of  the 
target  echo  and  the  false  target  are  less  than  the  resolution  of  the  radar,  interference  occurs. 
For  this  benchmark,  two  returns  are  asstuned  to  be  resolved  when  their  times  of  arrival  differ 


i  #. 


i  • 

I 

i 
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by  more  than  y.  Thus,  if 

C  T* 

the  target  and  RGPO  amplitudes  are  modified  according  to 

rl  =  n  +  [l  -  r{‘  cos  (2xe*) 

C  T« 

r{‘ = r{‘  +  fi  -  n  cos  (2t£*) 


c  T. 


(3.4) 


(3.5) 

(3.6) 


where  tk  =  —  H*  being  the  phase  of  the  false  target  signal,  and  c  is  the  speed  of 

light.  Note  that  the  discrete  coding  of  phase  has  been  largely  neglected  in  the  modeling  of 
the  interference  of  the  two  returns. 


The  RGPO  is  activated  by  a  radar  dwell  on  the  target.  Thus,  RGPO  reference  times  are 
dependent  on  the  dwell  times  of  the  radar.  For  the  benchmark,  RGPO  will  be  activated 
on  the  next  pulse  after  a  time  specified  in  the  trajectory  file.  The  reference  time  to  will  be 
chosen  as  the  dwell  time  minus  a  uniformly  distributed  random  number  between  0.1  and 
0.5  ps.  Thus,  ik  -  to  will  be  uniformly  distributed  between  0.1  and  0.5  ps.  The  amplification 
factor  7]  is  chosen  as  a  uniformly  distributed  number  between  1  and  4.  The  rates  for  both 
modes  have  been  selected  to  yield  a  2  km  displacement  in  range  of  the  false  target  in  20  s. 
Thus,  the  pull-off  rates  are  Vpo  =  100  m/s  and  opc  =  10  m/s^. 

Cover  pulse  is  a  RGPO  technique  that  uses  echoes  without  the  waveform  coding  and 
attempts  to  cover  the  actual  target  return.  When  cover  pulse  is  used,  a  noncoherent  detection 
occurs.  Typically  the  range  resolution  of  the  noncoherent  detection  is  rather  bad.  For  this 
benchmark,  cover  pulse  is  not  considered. 
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CHAPTER  4 

TARGET  TRAJECTORIES 


The  tergeto  will  exhibit  RCS  fluctuatioiis  according  to  the  Swerling  3  type  model  and 
perf<mn  aa  much  as  7  g  of  lateral  acceleration  and  2  g  of  lonptudinal  acceleration.  Target 
range  can  vary  from  20  to  120  km,  while  the  target  elevation  angle  varies  from  T  to  80°. 
Since  only  one  radar  face  it  used,  the  bearing  of  the  target  will  be  confined  to  ±60°.  The 
average  RCSs  of  the  targets  ate  large  enou^  so  that  average  SNRs  of  18  dB  are  achievable 
with  the  highest  energy  waveform.  The  SOJ  remains  at  ranges  greater  than  150  km  and 
performs  less  than  2  g  of  acceleration.  The  RGPO  can  be  employed  only  when  the  angle 
between  the  target  heading  vector  and  radar  range  vector  is  within  180°  ±  60°.  While  six 
target  trajectories  are  specified  in  this  section,  the  tracking  algorithm  should  be  designed  to 
handle  targets  satisfying  these  general  specifications. 

The  first  target  trajectory  is  shown  in  Figure  4.1  and  represents  a  large  aircraft,  such 
as  a  military  cargo  aircraft.  FVom  an  initial  range  of  80  km,  the  target  flies  on  a  constant 
course  with  a  speed  of  290  m/s  at  an  altitude  of  1.26  km  for  the  first  minute.  The  aircraft 
then  performs  a  mild  2>g  turn  and  continues  on  the  new  course  for  a  period  of  30  s,  where  a 
3-g  turn  is  made,  and  the  lurcraft  flies  away  from  the  radar  to  a  final  range  of  70  km.  The 
ia  4  m^,  RGPO  starts  after  15  and  40  s,  and  the  SOJ  is  within  2  degrees  of  the  target 
Line-of-Sight  (LOS)  from  82  to  93  s. 

The  second  target  trajectory  is  shown  in  Figure  4.2  and  represents  a  trajectory  that 
would  be  expected  by  a  smaller,  more  maneuverable  aircraft,  sudi  as  a  Learjet  or  other 
similar  high-performance  commercial  aircraft.  Target  2  initializes  at  a  range  of  63  km,  a 
speed  of  305  m/s,  and  an  altitude  of  4.57  km.  The  target  performs  a  2.5'g  turn  through  90 
degrees  of  course  change.  After  the  turn,  the  target  then  descends  gradually  to  an  altitude 
of  approximately  3.05  km,  A  4-g  turn  rolling  out  to  straight  and  level  flight  is  performed  at 
a  constant  speed  of  305  m/s,  and  the  trajectoiy  ends  near  a  range  of  28  km.  The  o^vt  i> 
2  m^;  RGPO  starts  after  12,  50,  and  95  s;  and  the  SOJ  does  not  approach  the  target  LOS. 
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The  trejectoriet  of  Targets  3  and  4  are  shown  in  Figures  4.3  and  4.4  and  represent  medium 
bombers  flying  at  high  speeds  with  good  maneuverability.  Target  3  baa  an  initial  speed  of  457 
m/s  and  flies  straight  and  level  for  the  first  30  s.  A  4-g  turn  is  then  performed  through  a  45- 
degree  course  change.  Straight  and  level,  nonaccelerating  flight  is  continued  for  the  next  30  s. 
A  second  4*g  turn  tbro\igh  a  90-degree  course  change  to  straight  and  level  flight  is  performed 
while  the  aircraft  decelerates  to  a  speed  of  274  m/s.  The  is  1.5  m^;  RGPO  starts  after 
30,  100,  and  130  s;  and  the  SOJ  does  not  approach  the  target  LOS.  Target  4  maintains  a 
speed  of  251  m/s  and  an  altitude  of  2.29  km  for  the  first  30  s.  A  4-g  turn  is  performed 
through  a  course  change  of  45  degrees.  After  another  30  s,  a  O-g  turn  is  performed  as  the 
throttle  is  increased  to  full  afterburner.  The  aircraft  pitches  up  and  climbs  to  an  altitude  of 
4.57  km.  Following  the  climb,  straight  and  level,  nonaccelerating  flight  is  maintained  for  the 
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Figure  4.4  l^Ajectory  for  Target  4 
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completion  of  the  trajectory.  The  is  1.7  m^,  RGPO  starts  after  35  and  65  s,  and  the 
SOJ  is  within  2  degrees  of  the  target  LOS  from  5  to  12  s  and  118  to  122  s. 

Targets  5  and  6  are  shown  in  Figures  4.5  and  4.6  and  represent  fighter/attack  aircraft. 
Thrget  5  is  initialized  in  a  thrusting  acceleration  at  an  altitude  of  1.5  km.  After  a  period 
of  30  s,  a  5-g  turn  is  performed  while  maintaining  full  throttle.  This  turn  is  followed  20  s 
later  by  a  7-g  turn.  Following  the  second  turn,  straight  and  level,  nonaccelerating  flight  is 
performed  for  30  s,  upon  which  a  6-g  turn  is  pertormed  concurrently  with  a  pitch  up  and  a 
climb.  After  an  altitude  of  4.45  km  is  reached,  straight  and  level,  nonaccelerating  flight  is 
flown  for  the  completion  of  the  trajectory.  The  is  1.2  RGPO  starts  after  5,  25,  and 
52  s;  and  the  SOJ  is  within  2  degrees  of  the  target  LOS  from  75  to  83  s.  Target  6  starts  at 
a  speed  of  426  m/s  and  an  altitude  of  1.55  km.  Constant  speed  and  course  are  maintained 
for  a  period  of  30  s  upon  which  a  7>g  turn  is  performed.  The  new  course  is  maintained  for 
another  30  s.  A  6-g  turn  is  performed  while  the  throttle  is  reduced  and  the  aircraft  is  nosed 
over  in  order  to  decrease  altitude.  After  a  final  altitude  of  0.79  km  is  obtained  along  with  a 
time  span  of  30  s,  another  6-g  turn  and  full  throttle  is  commanded.  After  approximately  30 
8,  a  7-g  turn  is  performed;  upon  completion  of  the  turn,  striught  and  level,  nonaccelerating 
flight  is  maintained  for  the  completion  of  the  trajectory.  The  >>  1*^  RGPO  starts 
after  15,  75,  and  102  s;  and  the  SOJ  is  within  2  degrees  of  the  target  LOS  from  69  to  76  s 
and  98  to  126  s. 
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CHAPTER  8 

TRACIONG  ALGORITHM 


Figure  5.1  ehows  the  input  and  output  of  the  tracking  algorithm.  At  any  revisit  time, 
as  many  as  five  dwells  of  any  combination  of  waveforms  in  Table  2.2  can  be  requested  by 
the  tracking  algorithm.  The  (xdkctkm  of  dwells  requested  at  a  given  time  by  the  tracking 
algorithm  will  be  referred  to  as  a  dwell  set.  The  tracking  algorithm  routine  will  receive 
for  each  dwell  in  the  requested  set  the  >od  for  each  range  bin  in  the  range 

gate,  with  Ra  above  the  commanded  SNR  threshold  1^.  Each  quantity  is  provided  to  the 
tracking  algorif  om  in  a  matrix  where  columns  1  through  5  correspond  to  dwells  1  through 
5  of  the  requested  set,  respectively.  Also,  indicators  of  coherent/noncoherent  detections 
(not  used  in  this  benchmark  problem)  and  the  waveforms  used  in  the  previous  radar  dwells 
are  given  to  the  tracking  algorithm.  Since  the  first  detection  occurs  on  a  search  dwell,  the 
tracking  algorithm  will  be  required  to  initialise  the  track  from  a  measurement  of  range,  and 
the  bearing  and  elevation  angles  of  the  beam  during  the  detection  (i.e.,  no  monopulse  error 
correction  in  the  first  measurement).  The  time  for  the  next  measurement  and  the  pointing 
commands  in  range  fa,  bearing  iumI  elevation  e§  are  computed  by  the  tracking  algorithm. 
The  tracking  (and  radar  management)  algorithm  also  requests  the  type  of  dwell  (i.e.,  search, 
track,  or  passive)  and  selects  a  waveform  (i.e.,  pulse  length)  for  active  dwells  from  Table  2.2. 
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Figure  5.1  Input  and  Output  of  the  IVacking  Algoritlun 
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CHAPTER  6 

CRITERIA  FOR  EVALUATION  OF  ALGORITHMS 

For  evaluatioo  of  the  tracking  algorithms,  each  algorithm  shall  maintain  tracks  on  the 
given  trajectories  with  a  maximum  track  loss  of  4%  and  no  indicati<m  of  the  target  type  or 
number.  A  track  is  declared  lost  if  the  error  in  the  estimated  position  of  the  target  is  greater 
than  1  two-way  beamwidth  in  angle  or  1.5  range  gates  (about  2362  m)  in  range.  Since  the 
phased  array  radar  will  operate  in  different  environments,  the  primary  measure  of  tracking 
performance  will  be  a  weighted  sum  of  the  average  radar  energy  per  second  and  the  average 
radar  time  per  second.  The  average  energy  per  second  i>  the  sum  of  energy  of  the  track 
dweUs  requested  by  the  tracking  algorithm  divided  by  the  number  seconds  of  the  trajectory. 
The  radar  energy  for  a  given  radar  dwell  is  computed  by  multiplying  the  transmit  power 
Pt  by  the  effective  pulse  width  r,.  The  average  radar  time  per  second  is  the  sum  of 
radar  time  of  the  track  dweUs  requested  by  the  tracking  algorithm  divided  by  the  number  of 
seconds  of  the  trajectory.  Each  radar  dwell  requires  0.001  s  of  radar  time.  Note  that  7^ 
is  inversely  proportional  to  the  average  sample  period.  The  quantities  and  are 
computed  for  each  Monte  Carlo  experiment.  Two  cost  functions  for  each  target  are  defined 
as 

Cl  =  (6.1) 

C2  =  +  10*7.,,  (6.2) 

where  Ci  corresponds  to  a  period  of  operation  when  radar  energy  is  critical,  and  C3  corre¬ 
sponds  to  a  period  of  operation  when  radar  time  is  critical.  Note  that  the  objective  is  to 
minimize  one  of  the  two  cost  functions  upon  request.  However,  Ci  and  C2  should  be  com¬ 
puted  for  each  target  and  presented  as  the  results  may  indicate  the  particular  benefits  of  a 
given  algorithm.  For  a  final  assessment  of  algorithm  performance,  each  cost  in  (6.1 )  and  (6.2) 
are  to  be  averaged  over  all  the  targets,  with  the  costs  from  Target  1  being  counted  twice. 
The  algorithm  that  provides  the  minimum  of  both  costs  by  a  simple  change  of  the  design 
parameters  will  be  considered  the  superior  algorithm.  Note  that  the  energy  and  radar  time 
during  initialization  should  be  included  because  initialization  is  a  portion  of  the  tracking. 
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Table  6.1  Summary  of  Results 


TWget 

No. 

Tim« 

Len 

(•) 

Max 

Accel 

(m/s*) 

Man 

Den 

(%) 

Sample 

Period 

.  .  _(•)  .  . 

RMSE 

Poa 

_  (m)  _ 

RMSE 

Vel 

(m/a) 

c, 

C2 

Track 

Lost 

.  .(%) 

1 

165 

31 

29 

2 

150 

39 

34 

3 

145 

42 

36 

i 

185 

58 

46 

5 

182 

68 

70 

6 

188 

70 

59 

Ave' 

*  lliii  average  is  computed  with  results  of  Target  1  being  counted  twice. 

A  secondary  measure  of  algorithm  performance  is  the  computer  resources  required.  Each 
investigator  will  present  an  assessment  of  the  average  number  of  floating  point  operations 
per  second  that  their  algorithm  will  require  in  final  implementation.  If  the  algorithm  has 
computational  requirements  that  vary  over  time,  an  assessment  of  the  peak  number  of  float 
point  operations  per  second  should  also  be  presented.  The  average  sample  periods  and  Root- 
Mean-Square  Errors  (RMSEs)  in  the  position,  velocity,  and  acceleration  estimates  should  be 
plotted  for  illustration  of  algorithm  performance.  The  results  should  be  reported  in  a  table 
as  shown  in  Table  6.1,  where  the  average  of  the  results  are  to  be  reported  in  the  last  row. 
The  fourth  coltunn  indicates  the  maneuver  density  which  was  computed  as  the  percent  of 
the  total  time  that  the  target  acceleration  exceeds  5  m/s^.  The  fifth  column  is  the  average 
sample  period,  while  the  sixth  and  seventh  columns  are  the  RMSE  in  position  and  velocity, 
respectively.  The  eighth  and  ninth  columns  are  the  costs  of  (6.1)  and  (6.2),  respectively.  The 
tenth  column  is  the  percentage  of  tracks  that  are  lost  during  the  Monte  Carlo  simulation. 
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CHAPTER  7 

CONCLUDING  REMARKS 


Whea  devdi^ing  »  sdutioii  to  this  baichro«rk  problem,  <me  common  tracking  algorithm 
should  be  designed  and  implemented  for  all  six  target  trajectories.  The  tracking  algorithm 
should  be  designed  to  address  the  general  specifications  of  the  targets  and  evaluated  with  the 
six  trajectories.  While  the  tracking  algorithm  can  adapt  to  the  trajectory,  the  adaptation 
should  be  automatic  in  that  no  user-defined  inputs  denote  a  target  type  or  specific  number. 
After  a  table  similar  to  Table  6.1  is  completed  to  summarize  the  results,  the  strengths  of  the 
solution  should  be  discussed  along  with  the  weaknesses.  Note  that  while  the  best  algorithm 
is  the  one  that  minimizes  the  two  cost  functions  of  Section  6,  the  computational  resources 
required  to  achieve  a  given  performance  are  very  important  when  selecting  a  tracking  method 
for  a  radar  system.  Thus,  solutions  with  various  computational  requirements  are  of  interest. 
FWthermore,  this  benchmark  problem  can  be  used  to  illustrate  the  application  of  different 
tracking  methods.  The  computer  progranu  that  are  enclosed  in  electronic  format  allow  the 
effects  of  the  RGPO  and/or  the  SOJ  to  be  easily  removed  from  the  simulation  program  by 
a  simple  change  of  the  input  file,  as  discussed  in  Appendix  B. 

In  the  development  of  this  benchmark  problem,  a  number  of  simplifying  assumptions 
were  made  in  order  to  limit  the  scope  of  the  problem  and  the  complexity  of  the  simulation 
program.  First,  the  targets  were  modeled  as  point  targets  with  RCS  fluctuations  that  were 
independent  of  the  aspect  angle  of  the  target  with  respect  to  the  radar.  In  an  actual  radar 
tracking  system,  the  returns  from  the  targets  can  occupy  multiple  range  bins,  glint  errors  are 
common  for  targets  at  the  closer  ranges,  and  the  RCSs  of  targets  can  change  rapidly  when 
targets  maneuver.  Second,  the  effects  of  closely  spaced  targets  have  not  been  considered.  If 
two  targets  are  separated  by  less  than  a  beamwidth  in  angle  and  their  ranges  are  not  fully 
resolved,  the  returns  from  the  two  targets  will  interfere.  In  a  monopulse  tracking  system,  this 
interference  can  be  catastrophic  to  the  tracking.  Third,  sea-surface  induced  multipath  has 
been  neglected  by  considering  only  targets  with  altitudes  above  the  region  where  multipath 
reflections  corrupt  the  monopulse  processing.  Fourth,  the  power  of  the  SOJ  has  been  limited 
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80  that  it  can  be  defeated  with  the  higher  energy  waveforms.  In  an  actual  system,  the  SOJ 
power  could  require  the  averaging  of  the  returns  from  many  dwells  or  coasting  the  track, 
while  both  the  target  and  the  SOJ  are  in  the  mainbeam.  Also,  if  the  SOJ  power  is  sufficiently 
high,  sidelobe  jamming  can  become  a  serious  issue.  Fifth,  the  potential  for  a  jammer  onboard 
the  target  has  been  ignored.  Sixth,  the  waveform  types  have  been  limited  to  fixed  waveforms 
with  discrete  codes.  Pulse  Doppler  and  linear  Frequency  Modulated  (FM)  waveforms  are 
commonly  used  radar  waveforms,  while  adapting  the  discrete  coding  of  the  waveforms  may 
prove  beneficial  in  the  presence  of  unresolved  targets,  multipath,  or  clutter.  Seventh,  the 
effects  of  background  clutter  have  been  neglected.  Eighth,  track  initiation  in  a  cluttered 
environment  has  not  been  considered.  Many  of  this  issues  are  open  research  problems  to  be 
considered  in  the  future. 
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APPENDIX  A 

DERIVATION  OF  MONOPULSE  PROCESSING  EQUATIONS 
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While  the  problem  of  tracking  multiple  targets  has  been  studied  extensively  in  recent 
years,  the  issue  of  finite  sensor  resolution  has  been  completely  ignored  in  almost  all  of  the 
studies  [A-1].  Typically,  the  targets  are  assumed  to  be  detected  with  a  given  probability 
of  detection  in  the  presence  of  false  alarms  and  clutter,  and  the  target  measurements  are 
modeled  as  the  true  values  plus  independent  Gaussian  errors  [A-2,A-3].  However,  when  the 
angular  position  of  a  target  in  the  presence  of  a  jammer  is  measured  using  a  monopulse  radar, 
the  Direction-Of-Arrival  (DOA)  information  for  the  target  and  jammer  are  merged  into  one 
measurement.  DOA  estimation  for  a  target  in  the  presence  of  a  Gaussian  noise  Standoff 
Jammer  (SOJ)  is  considered  in  this  Appendix.  First,  some  background  on  monopulse  systems 
and  related  processing  is  given  along  with  definitions  of  the  associated  notation.  Second,  the 
conditional  probability  density  function  (pdf)  and  related  statistical  properties  of  monopulse 
measurements  of  a  target  in  the  presence  of  a  Gaussian  SOJ  are  considered.  Third,  DOA 
estimation  for  an  SOJ  in  the  absence  of  a  target  is  developed.  Fourth,  the  DOA  estimation 
is  developed  for  a  target  in  the  absence  of  an  SOJ.  Fifth,  DOA  estimation  for  a  target  in  the 
presence  of  a  SOJ  is  developed  utilizing  the  estimates  of  Jammer-to-Noise  Ratio  (JNR)  and 
DOA  for  the  SOJ  that  can  be  formed  with  the  measurements  in  the  range  gate  that  do  not 
include  the  target. 

Background  and  Definitions 

i' 

In  a  typical  monopulse  radar  system,  the  outputs  of  the  receivers  are  matched  filtered, 
and  the  in-phase  and  quadrature  portions  of  the  sum  and  difference  signals  for  the  merged 
measurements  of  a  target  and  SOJ  can  be  expressed  as 

sj  =  aj  cos<l>j  +  at  cos<i>t  +  nsj 
SQ  —  Oj  +  at  sin^i  +  nsQ 

dj  =  ajr\j  cos<f>j  +  attjt  cos^t  +  nn 
dq  =  ajTfj  sin^j  +  atr}t  sin^i  +  njQ 

where 


Qj  =  ^jGi:{9j) 

(A.5) 

at  =  VKAtG^{9t)pQ 

(AS) 

=  DOA  of  the  SOJ 

(A.7) 

GA(0i\ 

Vt  =  77-7/1  V  =  DOA  of  the  target 

(A.8) 

K  =  proportional  to  the  transmitted  power 
Aj  —  voltage  amplitude  of  the  SOJ 
At  =  voltage  amplitude  of  the  target  backscatter 
=  sum  channel  antenna  gain  at  the  angle  9 
~  difference  channel  antenna  gain  at  the  angle  9 


(A.1) 

(A.2) 

(A.3) 

(AA) 


A-3 
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Bj  s:  off-boresight  angle  of  the  SOJ 
Bi  ae  off-horesight  angle  of  the  target 
Pd  =  matched  filter  gain 
=  phaae  of  the  SOJ  signal 
=  phase  of  the  target  return  signal 
"5/  ^ (0,  <t|)  -v  fi (0,  <r|) 

nu  iV(0,  (f\)  n^Q  -w  iV(0,  ol) 

The  ngi,nsQ,nd/,  and  n^Q  are  assumed  independent.  Since  the  SOJ  transmits  Gaussian 
noise,  the  amplitude  of  the  SOJ  Aj  will  be  Rayleigh  with  parameter  Ajo,  while  the  amplitude 
of  the  target  At  will  be  assumed  fixed  since  only  ooe  observation  of  the  target  amplitude  is 
assumed  to  be  available.  Then  the  amplitude  of  the  SOJ  signal  at  the  radar  is  given  by 

/(ojlojo)  =  ^  exp  ,  oijo  =  AjoGz(Bj)  (>1.9) 

20^0. 


Letting  A  and  0  denote  the  measured  amplitude  and  phaae  of  the  sum-signal  gives 

s/  s  A  cos  V*  SQ  =  A  sin  V*  (>1.10) 

where  -r  <  il>  <  K.  The  pdf  of  A  is  found  by  applying  the  transformation  of  random 

variables  in  (A.  10)  to  the  pdf  of  s/  and  tq.  Since  aj  is  Rayleigh  distributed  and 

uniformly  distributed  on  (*-x,x],  s/  and  s<}  are  independent,  Gaussian  random  variables  for 
given  values  of  at  and  Thus, 

£(s/|0|  ss  oi  cos^  (A.ll) 

£(s<]>|6]  =  at  sin^  (>1.1^) 

where  0  denotes  the  parameter  set  The  variances  are  given  by 

VAR(s/|0]=  VARIsq|0]  =  oJo  +  <t|  (A.13) 

Thus,  the  pdf  of  A  is  given  by 

~ 

where  /o(‘)  is  the  zero-order  modified  Bessel  function  of  the  first  kind. 


In  order  to  write  the  results  of  later  sections  in  terms  of  SNR  and  JNR,  let 


(A.15) 


A-4 
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where  J  rq>reseiits  the  JNR  of  the  SOJ,  I?  represeots  the  SNR  of  the  t&rget,  and  Rp  represents 
the  observ^  SNR. 


The  DOA  measurements  of  a  monopulse  radar  are  typically  taken  as  the  in-phase 
monopulse  ratio  (A-4].  With  s  =  s/  -f  jsQ  and  d  —  di  +  jdQ,  the  in-phase  and  quadra¬ 
ture  monopulse  ratios  (i.e.,  measurements)  are  given  by 


(Ai6) 

(A17) 


Distribution  of  Monopulse  Ratio 


Since  ay  is  Rayleigh  distributed  and  is  uniformly  distributed  on  (-r,  t],  the  sum  and 
difference  signals  are  Gaussian  distributed  given  at  and  The  pdf  of  the  sum  and  difference 
signals  is  defined  by  the  mean  X  =  E[X\0\  and  covariance  P  =  E[{X  -  ^)(A'  — 
whereA‘=[s/  dj  »q  Then 


Of  coe^ 

Pll 

P12 

0 

0  ■ 

at  sin^ 

P  — 

P12 

P22 

0 

0 

0|i7<  coi^ 

a  • 

0 

0 

Pll 

P12 

attjt  sin^i, 

0 

0 

P12 

P22. 

(A.18) 


where 


Pll  =  ojo  +  <^S 


Pl2  * 


Pn  = 

The  inverse  and  determinant  of  P  are  given  by 


(A.19) 

(A.20) 

(A.21) 


P-^  = 


(piiP22  -  pja) 
l^i  =  (PllPM  -  ^2)* 


P22 

“PI2 

0 

0 


-Pl2 

Pll 

0 

0 


0 

0 

P22 

“P12 


0 

0 

-’P12 

Pll 


(A.22) 

(A.23) 


Then 


X  exp 


{ ■  fe(‘"  ■  ■  “■ 
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(^.24) 


(>4.25) 

(4.26) 


+  ((d/  -  atm  cw^)*  +  (dQ  -  Oifft  sindf)*) 

2pia  / 

-  —  I  («/  -  co8^)(d/  -  iivH  coe^) 

Pii  ' 

+  (<Q  -  Qti  8invt)(dQ  -  atTii  rin^))  j  | 

Applying  the  transfonnntion  of  random  variables  in  (A.IO)  and 

dj  s  sjyi  -  aqyq  =  y/A  cos^  ~  yqS.  sind* 
dq  =  sjyq  +  $qyi  =  y^A  coed'  +  y/A  8ind» 

gives 

. 

-  2— ^A^y/  “  OfAy/  cos  e  <(-  oiAyo  sin  e  -  oiAiii  cos  e  +  ot(»7i) 

Pii  ' 

+  (A^(y?  +  Pq)  -  2y;iha|A  cos  c  +  2yqr)tat\  sin  e  +  j  |  (4.27) 

where  A  >  0,  -oo  <yi  <  +oo,  -oo  <yq  <  +oo,  and  -jt  -  d«  <£<*■-  di  with  £  -.  d*  ~ 
Then 


—  f  +  Oj  —  2a|A  COB 
Pll  '>  ' 


X  exp 

+ 

+ 


2(PiiP33  _ 

+  2~7«ai(A  COB  e  -  a<)  +  a?»?t  j  | 


Pll 

A» 

4*’*(PiiP32  -Pia) 


X  exp 

+ 

+ 


"Ox  "”')]} 


A-6 
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Using  (A.  14)  to  condition  the  pdf  of  (A.28}  on  A  gives 


X  expj 
X  expJ 


If  a|  s  0,  integrating  (A.29)  with  respect  to  e  gives 

ti  lA  a  m  A^pii  r  A^pii[y^  +  (yj-pi2pf,*)^]^ 

VUVQ  ,  ,o«=  “  2»(piipw  -;^j)  2(^iP23-Pi2)  ^ 

=  /(w|A,0)/(yo|A,e)  (A.30) 


where 


/(»,|A,e).^(^,EU^)  +  (4.31) 

/(v,|A,e)  =  Al(o,  eu^)  .  Al(o.^[g  +  7^,?])  (A.32) 


and  y)  denotes  the  Gaussian  distribution  with  mean  x  and  variance  y.  Thus,  for  the 
SOJ  in  the  absence  of  the  target,  yj  and  yq  are  conditionally  independent,  Gaussian  random 
variables  given  A. 

If  both  targets  are  on  the  antenna  boresight  (i.e.,  tfjssrjt^  0),  (A.29)  becomes 
fiVi, yg|A, e, s  17,  =  0)  =  /(y/|A, 0, 17,  =  17,=  0)/(y<j|A, 0, 17,  =  t7«  =  0)  (A.33) 

where 

A. 


/(w|A,0,»7,  =  i7i  =  O)  =  Ar(o,^) 
/(yQ|A,0,i7>  =  «7<  =  0)  S=  iv(o, 


(A.34) 

(A,35) 


Thus,  yj  and  yq  are  conditionally  independent  Gaussian  given  A  when  rfj  =  171  s  0.  However, 
in  general  (i.e.,  tjj  ^  0  or  171  7^  0  and  Of  ^  0)  /(WiVqIAiG)  has  a  difficult  form  for  computing 
the  moments  directly.  Thus,  the  joint  characteristic  function  [A-5]  is  used  and  given  by 

/•  r<»  roo 

/  /  exp[7u>iy/  ^  jW3yg]/(e,  y/,  y<j|A,  0)  dyj  yq  de 

't  J^OO  J-^OO 


A^Pij 


-Pi3)/o(o»ApJ 


— rr  /  dc  exp  — cos  c 

Ipii  -I 
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junvi  - 


Pn 


2(PhP32 
PliA* 


“■)))<" 


1 


2T/o(o|Apj7)  \  2piiA2 

-(w?+w^)  +  >~u;i| 

Pil  J 


{- 


(PiiPaa  ~  A) 


(w?  +w|) 


pii  J 


1 


/o(aiApn  ) 


rn-exp 


{- 


(piiPm  ~  Pia)/  2 
2piiAi 


X  /o  (x^ ~  (PiaPn  “  '^i))  “  “'a  (piaPn  “  ’?«)  ) 

The  marginftl  characteristic  functioos  of  j^/  and  are  then  given  by 

*  /  N  nx  1  ...  (  (P11P22-P12)  2  .  Piz  1 

*«(“■) = “  skxs?)  “n — ■''^“‘1 

«  A)(j(aW  -  j"iO>iW  -  'll))) 

X  lo  ^  (A*Pit^  -  (^l(pt2Pn  -  Vt)^ 


(A.37) 


(A.38) 


which  shows  that  y/  and  pQ  are  not  Gaussian  ruadoro  variables  if  a|  ^  0  and  rjt  ^  PiaPn^- 
However,  since  is  real  for  all  u^s,  the  pdf  of  yg  is  symmetrical  about  zero,  which 

gives  f^[yg|A,^]  =  0.  Furthermore,  since  for  Of^O  and  ijt  ^  PnPn^ 


(.4.39) 


yj  and  yg  are  not  in  general  independent  random  variables.  Setting  S]  =  ;u>i  and  »2  —  ji^ 
in  (A.36)  gives  the  joint  moment  generating  function  as 


M{»i,S2)  = 


1 


/o(a<Apf,*) 


PTT^'P 


{ 


(PiiPn  -j^)  1,2  ^  .2\ 


2^ 


+  4)  + 


r*'} 

Pii  ) 


X  “  *i(Pi3Pn  -  m)f  +  *2(PiJPn  ”  (>4.40) 

The  form  of  (A.40)  suggests  the  use  of  joint  cumuiants  for  the  computation  of  moments  of 
y;  and  yg.  The  joint  cumuiants  [A-d]  are  given  by 


.  I 


(>4.41) 
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where 


(PiiP23  —  Pia)/-!  .  -a\  .  Pi2 


*{sus2)  =  =  -In  [/o(a.Apr,‘)]  +  + 


— 

Pn 


*n  (a^Pu*  -  {PiiPu  ~  ’7«)«i)  +  (PiaPn  ”  (^-42) 


The  iiist'Order  moments  or  expected  values  of  yj  and  yq  are  given  by 


^  2 
^(y/|A,*l »  Aio  «  ^3  +<r|^  ~  y^A|o(a<A(o^o  +  f  ~ 

£lyQ|A,*)  =  Aoi=0  (AM) 


where 


(A.45) 


with  /i(-)  denoting  the  first-order  modified  Bessel  function  of  the  first  kind.  The  second 
central  moments  of  yj  and  yq  are  given  by 


^2 -a 


VARl„|A,e|  =  ^[^  +  '  ’”) 

X  I  -  /f|Q(o|A(<rjo  +  <^s)"*)  “  ^•~^^-]^~~A|o(Q<A(a^o  +  ‘^5)"')j  j  (A.46) 

vA%,|A,e|  =  g[^+jri^ 


cov(i,M,,|A,e|  =  o 


(a*o  +  ~  A|o(®‘A(o*o  +  <^|)”’) 


(A.47) 

(A.48) 


where  VAR|yf  j-]  denotes  the  variance  of  yj,  and  COV(y/,y^|’)  denotes  the  covariance  of  yj 
and  yq.  While  y;  and  yq  are  not  independent,  (A.48)  indicates  that  the  random  variables 
are  conditionally  uncorrelated.  Thus,  any  statistical  processing  of  y/  and  yq  that  is  based 
on  the  first  two  moments  can  be  accomplished  independently.  Farther  analysis  that  includes 
the  third-  and  fourth-order  marginal  cumulants  of  y/  and  yq,  and  the  kurtosis  as  conducted 
in  [A-fi]  for  the  case  of  a  single  fixed-amplitude  target  will  indicate  that  y/  and  yq  can  be 
closely  approjdmated  as  independent,  Gaussian  random  variables. 

DOA  Estimation  for  80J  in  the  Absence  of  a  Ikrget 

Since  the  energy  of  the  SOJ  enters  into  all  of  the  bins  of  the  range  gate,  all  of  the  bins 
can  be  used  to  estimate  the  DOA  of  the  SOJ.  Setting  Oi  =  0  in  (A. 14)  and  applying  the 
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transform&tion  of  (A. 15)  in  (A.14)  gives  the  distribution  of  the  observed  Signal-to-Noise 
Ratio  (SNR)  of  each  range  bin  for  measurements  of  the  SOJ  as 

/(»*l^)=  »o>0  (A49) 

For  N  samples  or  range  bins,  the  Maximum  Likelihood  (ML)  estimate  of  J  is  given  by 

J  =  (^-50) 


where 


(A51) 


with  SSoi  denoting  the  observed  SNR  for  bin  t.  The  estimate  J  is  unbiased  with  variance 
given  by 

/  J  I  1 \Z 

(A.52) 


(J  + 1) 

VAR[J|J1  = 


Since  y/,  and  yQi  are  Gaussian  distributed,  the  ML  estimate  of  y;,  the  mean  of  y/,  is 
given  for  N  range  bins  by 

is]  isl 

where  Sfoi  Vli  denote  the  observed  SNR  and  in-phase  monopulse  ratio  for  bin  i.  Thus, 
the  estimate  y/  is  a  “power”  weighted  sum  of  the  N  monopulse  ratios.  Since  each  y/,-  is 
a  Gaussian  random  variable,  y/  is  the  minimum  variance  estimate  of  yj  and  a  Gaussian 
random  variable,  with  variance  given  by 

N  J 


where 


(A.55) 


Using  the  Method  of  Moments  [A-7]  estimation  and  yj  to  estimate  i/;  gives 

7y=[l-hi]y/,  J>3dB  (A.56) 

where  JNR  has  been  replaced  by  its  ML  estimate  since  the  JNR  is  typically  unknown.  Note 
that  J  =  J  is  a  appropriate  for  N  >  8. 


Using  the  mean  of  y;  in  (A.31)  to  form  an  estimate  of  p  gives  an  estimate  of  the  variance 

N 

c 

i-1 


(A.57) 
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where 

^  ^  >  3  dB  (>1.58) 

Using  the  resulting  p,  an  estimate  of  the  variance  the  monopulse  ratios  of  bin  t,  can  be 
computed  using  the  N  bins  according 


x2  _  P 


(A.59) 


Using  (A.56)  and  (A.57)  gives  an  estimate  of  the  variance  associated  with  ifj  as 

+ j]  ^?/ =  ^  +  ^1  ^>3dB 


(A.60) 


DOA  Estimation  for  Ikrget  in  the  Absence  of  an  SOJ 

When  a  single  pulse  is  used  to  form  a  monopulse  ratio  for  DOA  estimation,  the  target 
amplitude  can  be  treated  as  fixed  and  the  phase  is  uniformly  distributed  in  the  interval 
(—X ,  x].  The  target  amplitude  is  observed  through  the  measured  amplitude  of  the  sum^signal, 
where  the  pdf  of  A  is  found  by  setting  aj  =  0  in  (A.  14).  Applying  the  transformations  of 
(A.  15)  in  the  pdf  gives 

/(»o|R)  =  /o(2v^)  exp{  -  («,  +  »)},  »„  >  0  (A.61) 

Then  the  ML  estimate  of  Hi  based  on  a  single  pulse  satisfies 


(A.62) 


where  /i|o(i)  is  defined  in  (A.45). 


The  conditional  pdf  of  the  yi,yQ,  and  e  is  found  by  setting  pis  =  0  in  (A.29)  to  obtain 

..  A*  exp{o,A<Tj2  cos  e}  ^  A*  /  q, 

/U,W.V«|A,e)  =  «cp{  -  cos  e)  } 


1  *’/ 

f  Of  .  \ 

^VQ  -  jm  e ) 

(A.63) 


Since  the  location  of  the  pdf  of  yq  is  independent  of  r/t  (j.e.,  £^(pq|A,©]  =  0),  a  single 
observation  of  yq  provides  essentially  no  information  concerning  the  value  of  rff.  Note  that 
while  (A.47)  shows  that  the  variance  of  yq  depends  on  »7j,  single-pulse  estimation  of  the 
variance  is  not  considered  viable.  Noting  that  the  ML  value  of  yq  is  zero  and  using  PQ  =  0 


A-11 


NSWCDD/TR-96/10 


in  (A.63)  gives  e  s  0.  Thus,  the  ML  estimate  of  rjt  for  a  known  target  amplitude  at  is  given 

(/1.64) 

Comparing  (A.64)  and  (A.43)  with  ayo  =  0  clearly  indicates  that  (A.64)  gives  a  biased 
estimate  of  171.  While  an  unbiased  estimate  of  r/t  can  be  computed  with  (A.43),  the  estimate 
will  not  be  ML.  Using  ^  and  a^o  =  0  in  (A.46)  and  (A.47)  gives  an  estimate  of  the  variance 
of  yj  and  yq  as 


1 

2R, 

A 

A2 


if  +  ^Vl^e 


/T-  /t2 

«2  _  ,  .,2  j  a^-2\ _ 


QfA 


2»o 


^;^2,|0(2v/^)]] 


(A.65) 

(A.66) 


Using  (A.65)  and  (A.66)  gives  an  estimate  of  the  variance  of  for  a  known  target  amplitude 
as  =  RoSi'^dJ^,  which  is  required  for  tracking  with  a  Kalman  filter  or 

Interacting  Multiple  Model  (IMM)  algorithm  [A-2]. 


When  the  target  amplitude  is  unknown  to  the  signal  processor,  a  ML  estimate  of  the  target 
amplitude  will  be  used  along  with  the  monopulse  ratio  to  estimate  the  DOA.  Since  the  target 
amplitude  is  observed  through  the  sum-signal,  the  estimation  of  the  target  amplitude  will 
be  decoupled  from  the  monopulse  processing.  Given  an  observed  SNR  IRo<  (A.62)  gives  the 
ML  estimate  of  Note  that  S?o  is  a  biased  observation  of  91.  An  approximate  ML  estimate 
is  given  by 


where  the  maximum  percent  error  is  5  percent  at  Slg  =  3  dB.  The  limitation  of  ^  3  dD 
is  not  a  particular  problem  since  the  detection  threshold  is  often  greater  than  7  dB  and  a 
single  observation  with  a  smaller  amplitude  provides  essentially  no  error  reduction  relative 
to  the  beamwidth  of  the  sum  pattern.  Thus,  using  %  for  91  in  (A.64)  gives  an  approximately 
ML  estimate  of  i/i  for  a  monopulse  measurement  of  target  with  unknown  amplitude  as 


91. 


91o  +  0.5 


91.  >  3  dB 


(A.67) 


m  = 


4 


I  _ ^ 


y/  = 


1  + 


291. 


VI 


91.  >  3  dB 


(A.68) 


While  (A.68)  provides  a  biased  estimate  of  tft,  the  bias  will  be  less  than  that  provided  by 
standard  monopulse  ratio.  A  nearly  unbiased  estimate  can  be  computed  by  using  (A.62)  and 
(A.67)  in  (A.43)  with  o^q  =  0,  but  the  estimate  will  not  be  ML.  Under  the  ML  conditions 
of  (A.62) 

^  ^  1  -  (^-69) 
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Using  (A.69)  in  (A.65)  and  (A.62)  in  (A.66)  with  li »  ^  gives  estimates  of  the  variances  of 

yi  and  yq  as 


A2  ^  2A« 


.JlL 


38,  >  3  dB 
38,  >  3  dB 


(A.70) 

(A.71) 


Note  that  differs  significantly  from  that  reported  in  [A-4],  which  is  restricted  to  high 
SNRs.  Monte  Carlo  simulation  studies  have  been  conducted  for  various  values  of  38  and  rjt 
to  assess  the  validity  of  (A.70)  for  estimating  the  variances  of  the  corresponding  monopulse 
ratios.  For  each  experiment  the  amplitude  was  restricted  so  that  38,  >  3  dB  and  <rj  =  <r|, 
and  yi  and  yg  were  normalized  by  the  corresponding  standard  deviation  of  (A.70)  before 
computing  the  sample  standard  deviation.  The  simulation  results  indicc;.ed  that  the  variance 
estimate  of  yj  is  v4d  for  3  <  38  <  12  dB  with  |i;(|  <  2  and  for  38  >  15  dB  with  i>7(|  <  3. 
Note  that  the  variance  given  in  (A-4]  was  found  valid  for  38  >  12  dB  and  |>7(|  <  0.25.  The 
results  for  yq  indicated  that  is  valid  for  3  <  38  <  9  dB  writh  |>7|  <  1  and  38  >  12  dB  with 
l»ll  <  3. 


Using  (A.68)  and  (A.70)  gives  an  estimate  of  the  variance  of  m 


= 


2 

<’3  .  t>M 

1 

[i-.  M 

2 

,  Vj 

1  4. 

1 

“  2», 

^■^  288, 

[;j+4s;] 

88,>3dB  (/1.72) 


which  is  required  for  tracking  a  target  with  an  unknown  amplitude.  Note  that  this  result 
differs  significantly  from  (A-4]  at  low  to  moderate  SNRs. 


DOA  Estimation  for  Target  in  the  Presence  of  a  SOJ 

When  estimating  the  DOA  of  a  target  in  the  presence  of  a  SOJ,  the  bins  within  the 
range  gate  that  do  not  include  the  target  can  be  used  to  estimate  the  J  and  rfj  with  (A.50) 
and  (A. 56)  as  discussed  above.  The  estimates  of  J  and  t)j  will  then  be  used  with  (A.29)  to 
estimate  the  DOA  of  the  target  and  (A.46)  to  estimate  the  associated  variance 

Applying  the  transformations  of  (A.  15)  to  the  pdf  of  (A.14)  gives 

/(»«!»)  «  /o(2v^(J  +  l)-‘)  exp{  -  K  >  0  (A.73) 

Then  the  ML  estimate  of  38,  based  on  a  single  pulse  given  that  J  is  known,  satisfies 


(A.74) 


A-13 


NSWCDD/TR.96/10 


Given  an  observed  SNR  Ko«  (A. 74)  gives  the  ML  estimate  of  9?  given  that  J  is  known.  An 
approximate  ML  estimate  of  is  given  by 


(^.75) 


«  r  ««  1  «  «  r  « 


Also,  under  the  ML  conditions  of  (A.74), 


Considering  J  and  rjj  as  known  in  (A.29)  ^ves  a  pdf  for  computing  the  ML  estimate  of 
ijt.  Since  the  location  of  the  pdf  of  yq  is  independent  of  rjt  (i.e.,  £[yQ|A,6]  =  0),  a  single 
observation  of  yq  provides  essentially  no  information  concerning  the  value  of  rff.  Note  that 
while  (A.47)  shows  that  the  variance  of  yq  depends  on  tjt,  single-pulse  estimation  of  the 
variance  is  not  considered  viable.  Noting  that  the  ML  value  of  yq  is  zero  and  using  yg  =  0 
in  (A.29)  gives  e  =  0.  Thus,  the  ML  estimate  of  tjt  in  the  presence  of  a  SOJ  for  a  known 
target  amplitude  Of«  is  given  by 


J  ^  Ar  J  ] 


(.4.77) 


When  the  target  amplitude  is  unknown,  using  (A.75)  in  (A. 77)  gives 

Note  that  only  the  mean  of  the  monopulse  ratio  for  the  SOJ  appears  in  (A.78).  Using  M 
range  bins  that  do  not  include  the  target  return  to  estimate  J  smd  the  mean  of  the  monopulse 
ratio  for  the  SOJ  gives 


where 


=  ^  +  R<,>2J  +  2 


^  =  SiL^2i!(£i 

l^kml  ^ok 


1  " 


(A.79) 


(A.80) 


(A.81) 


with  yji  and  Roi  denoting  the  monopulse  ratio  and  observed  SNR  for  bin  t  that  does  not 
include  the  target  return,  and  y;  and  denoting  the  monopulse  ratio  and  observed  SNR 
for  the  range  bin  that  includes  the  target  detection. 
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If  J  and  are  treated  as  known  in  (A.78),  then  the  variance  of  can  be  related  to  the 
variance  of  yj,  which  is  given  by  (A.46).  Using  the  approximation  of  (A.76)  in  (A.4€)  and 
(A. 77)  in  (A.46)  and  treating  J  and  as  known  in  (A.78)  gives  an  estimate  of  the  variance 
of  as 


Using  M  range  bins  that  do  not  include  the  target  return  to  estimate  J  and  rjj  gives 

where  y}  and  J  are  given  by  (A.80)  ai:«l  i). 


Eq.  (A. 79)  provides  an  estimate  of  the  DOA  of  a  target  in  the  presence  of  an  SOJ  and 
(A.83)  provides  an  estimate  of  the  variance  associated  with  the  DOA  estimate.  However, 
note  that  rfj  and  J  have  been  used  as  known  quantities.  Then  a  sufficiently  large  number  of 
range  bins  should  be  used  to  estimate  J  and  tfj.  A  reasonable  criteria  for  M  is  MJ  >  18  dB. 
For  example,  an  SOJ  with  J  =  6  dB  would  require  about  16  range  bins. 
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A  simulation  testbed  for  the  Benchmark  Problem  has  been  developed  by  G.  A.  Watson. 
The  simulation  is  written  in  MATLAB  for  version  4.0.  A  copy  of  the  simulation  code  has 
been  provided  with  this  report  on  computer  diskette.  A  brief  explanation  of  each  function  of 
the  simulation  program  is  given  below. 

1.  abaJ.ad.B  calculates  average  absolute  tracking  error  for  an  angle-only  Alter. 

2.  avrg-iad.B  calculates  average  track  Alter  output  values. 

3.  b«B»dth2 .  ■  calculates  the  bearing  and  elevation  beamwidths  for  the  given  radar  model. 

4.  bnfel.n  clears  the  buffer  when  displaying  the  status  report. 

5.  clcspace.n  clears  the  screen  and  prints  blank  lines. 

6.  coui2a.n  gates  passive  dwell  measurements. 

7.  cou2n.B  gates  track  and  search  dwell  measurements. 

8.  coord  .a  converts  positions,  velocities,  and  accelerations  in  Cartesian  coordinates  to  the 
corresponding  spherical  coordinate  values. 

9.  corrolr  .B  constructs  a  correlated  measurement  error  covariance  matrix. 

10.  covplot.B  plots  the  diagonal  elements  of  the  error  covariance  for  a  Cartesian  track 
Alter. 

11.  CTa«s96.B  is  the  primary  function  for  testbed  simulation. 

12.  cvf iltd.B  is  the  tracking  routine. 

13.  fcv.B  creates  the  constant  velocity  state  transition  matrix. 

14.  findq.B  creates  the  process  noise  covariance  matrix. 

15.  findra.B  creates  an  uncorrelated  measurement  error  covariance  for  a  range,  bearing, 
and  elevation  measurement  vector. 

16.  gatBda.B  performs  the  measurement  update  for  an  angle  only  track  Alter. 

17.  gatBdf  .B  performs  the  measurement  update  for  a  track  Alter  using  track  and  search 
dwell  measurements. 

16.  gcv.B  creates  the  constant  velocity  process  error  matrix. 

19.  havala .  b  converts  the  predicted  state  estimate  from  Cartesian  to  spherical  coordinates. 
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20.  bprcalca.B  CAloilates  the  H  matrix  for  an  extended  Kalman  track  filter  for  three 
coordinates. 

21.  initaper.B  calculates  the  initial  state  estimate  using  a  least-squares  technique  for 
Cartesian  coordinates. 

22.  initlbe .  m  calculates  the  initial  state  estimate  using  a  least-squares  technique  for  angle- 
only  tracking. 

23.  IdcTaasO.ais  the  default  load-in  file  for  cvaes96. 

24.  Mkevect  .■  makes  a  vector  of  uniform  random  numbers. 

25.  ■alcv«ct2  .■  makes  a  vector  of  normal  random  numbers. 

26.  pointcT.n  determines  where  to  point  the  radar  for  the  next  dwell  using  a  constant 
velocity  model. 

27.  radaod4a.nis  the  radar  model. 

28.  randeoi.a  determines  random  intervals  of  RGPO. 

29.  mseplot.n plots  the  Root- Mean-Square- Error  (RMSE)  of  a  Cartesian  track  filter. 

30.  calcxilates  the  average  RMSE. 

31.  signaper  .■  calculates  the  initial  error  covariance  matrix  using  a  least-squares  technique 
for  Cartesian  coordinates. 

32.  signcobe .  ■  calculates  the  initial  error  covariance  matrix  using  a  least-squares  technique 
for  angle  only  tracking. 

33.  space. ■  prints  blank  lines. 

34.  stdvplot  .■  plots  the  standard  deviations  of  the  error  covariance  for  position,  velocity, 
and  acceleration  for  Cartesian  track  filters. 

35.  tiner.n  displays  how  long  a  simulation  will  last. 

36.  tlBeup.n  is  the  time  update  for  a  track  filter. 

37.  trakerr2  .B  determines  whether  the  estimation  error  is  above  some  specified  threshold. 

38.  voltpat2.B  computes  the  radar  voltage  pattern. 

39.  zconvert  .b  converts  the  radar  measurements  to  Cartesian  coordinates. 
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Six  target  trajectories  have  been  provided.  The  targets  are  target! ,  targ«t2,  tar¬ 
gets,  target4,  targets  and  targets  and  are  explained  in  Chapter  4.  There  is  also  a 
Standoff  Jammer  (SOJ)  trajectory  called  janer.  The  target  trajectories  are  generated 
with  a  20  Hz  data  rate  and  have  the  form; 

[xxxyyyzz  it  rcs^^  (F.l) 

where  t  is  the  trajectory  time,  rcs«M  the  average  radar  cross  section  of  the  target,  and 
/ecm  >>  ^  indicator  of  on-board  target  ECM  employed  during  flight.  The  SOJ  trajectory 
does  not  have  rcsave  hem-  '^he  state  estimates  from  the  target  track  Alter  have  the  form 

x  y  y  z  z  (B.2) 

and  the  state  estimates  from  the  SOJ  track  filter  have  the  form 

=  4  c  e  if  (B.3) 

The  target  and  SOJ  state  estimates  may  be  saved  in  the  form 

Xl'fl^lx  X  X  y  y  y  z  z  z  tf  (B.4) 

and 

Xjj  =  {4  4  4  c  e  e  if  (B.5) 

which  includes  acceleration  estimates  if  desired.  The  frinction  that  calculates  the  averages 
output  values  senses  the  size  of  the  estimate  state  vector  and  perfonns  the  calculations 
accordingly. 

The  function  to  run  the  testbed  simulation  is  cva«a96.  It  has  a  default  load-in  file  that 
explains  each  parameter  needed  for  the  track  filters.  Only  output  arguments  are  needed  to 
run  cvaeeOA  when  using  the  load-in  flle.  The  filter  parameters  are  set  in  the  load-in  file.  The 
tracking  algorithm  can  also  be  run  using  the  call-line  in  cvaa896 .  This  is  accomplished  by 
passing  in  all  the  required  parameters  specified  in  the  load-in  file  in  the  proper  location  in 
the  call-line.  For  ease  of  using  the  call-line,  the  parameters  are  passed  to  cvaes96  in  the  same 
manner  as  the  output  line  of  the  default  load-in  file.  The  advantage  of  using  the  call-line  is 
that  a  batch  job  can  be  run  by  varying  different  input  pvameters. 

There  are  several  outputs  from  cva«a96.  If  the  SOJ  is  not  employed,  there  will  be 
outputs  only  for  the  target.  When  the  SOJ  is  used,  outputs  for  both  the  target  and  SOJ  will 
be  given.  The  outputs  of  cvaeB99  are  given  below. 

1.  zyz.targ  is  the  target  trajectory. 

2.  zyz.jaa  is  the  SOJ  trajectory. 

3.  ave..f  ilt.targ  is  the  average  target  filtered  state  estimates. 
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4.  avaJilt.jaa  is  the  average  SOJ  filtered  state  estimates. 

5.  av«.xma«.targ  is  the  average  target  RMSE  with  the  first  row  being  position,  the  second 
row  velocity,  and  the  third  row  time.  If  the  state  vector  has  acceleration  estimates,  the 
third  row  is  acceleration  and  the  fourth  row  time 

6.  av«-abs.jas  is  the  average  SOJ  absolute  track  error  (ATE)  with  the  first  row  being 
bearing,  the  second  row  bearing  rate,  the  third  row  elevation,  the  fourth  row  elevation 
rate,  and  the  fifth  row  time.  If  accelerations  are  calculated,  the  third  row  is  bearing 
acceleration,  the  fourth  elevation,  the  fifth  row  elevation  rate,  the  sixth  row  elevation 
acceleration,  and  the  seventh  row  time. 

7.  ave^ov.targ  is  the  average  of  the  diagonal  elements  of  the  target  filtered  error 
covariance. 

8.  avajcov.jan  is  the  average  of  the  diagonal  elements  of  the  SOJ  filtered  error  covariance. 

9.  avaJt.targ  is  the  average  target  sample  period  with  the  first  row  being  the  sample 
period  and  the  second  row  time. 

10.  aveJt-jaa  is  the  average  SOJ  sample  period  with  the  first  row  being  the  sample  period 
and  the  second  row  time. 

11.  ct.targ  is  the  number  of  occasions  each  target  state  was  chosen  with  the  first  row  being 
the  frequency  and  the  second  row  time. 

12.  ct.jan  is  the  number  of  occasions  each  SOJ  state  was  chosen  with  the  first  row  being 
the  frequency  and  the  second  row  time. 

13.  ava^ter.targ  is  the  average  number  of  target  updates  per  run  for  all  maintained 
tracks. 

14.  ava-itar-jan  is  the  average  number  of  SOJ  updates  per  run  for  all  maintained  tracks. 

15.  vl.targ  is  a  matrix  specifying  when  and  why  target  tracks  were  declared  lost.  If  there 
are  lost  target  tracks,  the  first  row  is  the  Monte  Carlo  iteration,  the  second  specifies 
why  the  track  was  lost  (lack  of  measurements  =  1,  tracking  error  was  above  =  2),  and 
the  third  tow  is  the  declaration  time. 

16.  wl.jaa  is  a  matrix  specifying  when  and  why  SOJ  tracks  were  declared  lost.  If  there  are 
lost  SOJ  tracks,  the  first  row  is  the  Monte  Carlo  iteration,  the  second  specifies  why  the 
track  was  lost  (lack  of  meastuements  a  1,  tracking  error  was  above  =  2),  and  the  third 
row  is  the  declaration  time. 

17.  parc.targ  has  two  elements  with  the  first  being  the  percentage  of  target  measurements 
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detected  and  the  second  the  percentage  of  target  measurements  within  the  tracking 
g»te. 

18.  perc.jaa  has  two  elements  with  the  first  being  the  percentage  of  SOJ  measurements 
detected  and  the  second  the  percentage  of  SOJ  measurements  within  the  tracking  gate. 

19.  aTOjdv.targ  is  the  average  dwell  time,  look  energy,  and  power  for  each  target  dwell 
time. 

20.  p-crit  is  a  performance  criteria  vector  with  nine  elements  and  has  the  form:  1)  average 
target  sample  period,  2)  average  dwell  time  per  dwell,  3)  average  energy  per  dwell, 
4)  average  dwell  time  per  track  length,  5)  average  energy  per  track  length,  6)  average 
dwell  time  per  run,  7)  average  energy  per  run,  8)  average  position  error  per  dwell,  and 
9)  average  velocity  error  per  dwell. 

The  output  perc-targ  and  perc.jaa  are  outputs  specifically  designed  for  this  applica¬ 
tion  to  provide  further  performance  evaluation  of  the  algorithm.  You  may  wish  to  add  ones 
pertinent  to  your  track  algorithm.  When  adaptively  sampling  the  trajectory,  the  output 
matrices  will  be  very  large  because  the  same  sequence  of  measurements  will  not  be  chosen 
each  time.  As  a  result,  a  running  total  is  kept  of  how  many  times  each  trajectory  state  was 
sampled.  Since  the  output  matrices  are  large,  you  may  want  to  reduce  their  size  by  binning 
the  outputs  over  a  specified  time  interval,  such  as  each  0.5  seconds.  The  output  matrices  will 
contain  values  if  that  trajectory  state  was  chosen  at  least  one  time  for  a  maintuned  track. 
If  the  state  was  not  chosen,  there  will  be  no  values  in  the  output  matrices  corresponding  to 
that  time.  The  binning  is  accomplished  by  averaging  the  outputs  over  each  specified  time 
intervd  as  determined  by  the  time  values  of  the  output  matrix.  There  may  be  different  num¬ 
bers  of  outputs  in  each  interval  depending  on  whether  or  not  that  state  was  chosen  during 
the  simulation.  You  may  wish  to  average  all  the  values  in  the  time  interval  or  do  a  weighted 
average  based  on  the  frequency  each  state  in  the  interval  was  sampled.  This  will  allow  for 
the  trend  to  be  seen  without  having  to  plot  each  output  time. 

If  you  have  any  question  about  the  programs  or  how  to  run  them,  please  feel  free  to 
contact  Greg  Watson  anytime!!  Phone:  (540)  653-7378,  Email:  gawatsoOnswc.navy.mil. 
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